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Abstract. Wc here investigate the well-posedness of a networked integrate-and- 
fire model describing an infinite population of neurons which interact with one 
another through their common statistical distribution. The interaction is of the 
self-excitatory type as, at any time, the potential of a neuron increases when some 
of the others fire: precisely, the kick it receives is proportional to the instantaneous 
proportion of firing neurons at the same time. From a mathematical point of view, 
the coefficient of proportionality, denoted by a, is of great importance as the 
resulting system is known to blow-up for large values of a. In the current paper, 
we focus on the complementary regime and prove that existence and uniqueness 
hold for all time when a is small enough. The critical value for existence and 
uniqueness is made explicit in some popular examples found in the neuroscience 
literature and is also compared with numerical experiments. 

Keywords: McKean nonlinear diffusion process; renewal process; density esti- 
mates; integrate-and-fire model; network of neurons; neuroscience. 

1. Introduction 

The stochastic integrate-and-fire model of the membrane potential V across a 
neuron in the brain has received a huge amount of attention since its introduction 
(see [18] for a comprehensive review). The central idea is to model V by thresh- 
old dynamics, in which the potential is described by a simple linear (stochastic) 
differential equation up until it reaches a fixed threshold value Vf, at which point 
the neuron emits a 'spike'. Experimentally, at this point an action potential is ob- 
served, whereby the potential increases very rapidly to a peak (hyperpolarization 
phase) before decreasing quickly to a reset value (depolarization phase). It then 
relatively slowly increases once more to the resting potential (refractory period). 

Since spikes are stereotyped events, they are fully characterized by the times 
at which they occur. The integrate-and-fire model is part of a family of spiking 
neuron models which take advantage of this by modeling only the spiking times 
and disregarding the nature of the spike itself. Specifically, in the integrate-and-fire 
model we observe jumps in the action potential as the voltage is immediately reset to 
a value Vr whenever it reaches the threshold Vp, which is motivated by the fact that 
the time period during which the action potential is observed is very small. Despite 
its simplicity, the integrate-and-fire model has been able to predict the spiking times 
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of a neuron with a reasonable degree of accuracy [TUl [TT] , and as such provides the 
basis for a good first study. 

Many extensions of the basic integrate-and-fire model have been studied, including 
ones in which attempts are made to include noise and to describe the situation when 
many integrate-and-fire neurons are placed in a network and interact with each other. 
Indeed, in [T3] and |14j the following equation describing how the potential V; of the 
i-th neuron in a network of N behaves in time is proposed: 

j n j =-oo jjti 

for Vi(t) < Vp and where Vi(t) is immediately reset to Vr when it reaches Vp. Here 
Ii Xt {t) represents the external input current to the neuron, rji(t) is the noise (a white 
noise) which is importantly supposed to be independent from neuron to neuron, and 
the constants A, /3, a and a are chosen according to experimental data. Moreover, 
the interaction term is described in terms of t n , , which is the time of the n-th spike 
of neuron j, and the Dirac function 5. Precisely, it says that whenever one of the 
other neurons in the network spikes, the potential across neuron i receives a 'kick' 
of size ol/N. 

In the case of a large network i.e. when N is large, many authors approximate the 
interaction term in the above finite system by an instantaneous rate v(t), the so- 
called mean-firing rate (see for example [U [21 dH US])- However, in the neuroscience 
literature little attention is paid to how this convergence is achieved. Mathematically 
the mean-field limit as N — » oo must be taken, and it is not at all clear that one 
observes a propagation of chaos phenomenon. This is the subject of a forthcoming 
manuscript [5]. 

Just as important is the question of whether solutions exist for the resulting non- 
linear equation when the mean-field limit is taken, and if so what type of solutions 
are admissible. Again this is not treated rigorously in the neuroscience literature, 
and it is in this diction that we proceed here. Since the mathematical difficulties lie 
within the jump interaction term, we in fact isolate it by making two simplifications. 
Indeed, we suppose that there is no external input current (I^ xt (t) = 0), and that the 
interaction term is composed solely of the jump or reset part ((3 = 0). Although this 
is a big simplification from a neuroscience perspective, it still captures all the math- 
ematical complexity of the resulting mean-field equation, and standard approaches 
(see for example [20] ) would then allow extensions to more realistic models. 

Without loss of generality, we also take the firing threshold Vp = 1 and the reset 
value Vr = for notational simplicity. We finally choose to describe the problem 
in the language of stochastic processes, so that the precise nonlinear stochastic 
differential equation under study here becomes 

X t = X + [ b(X s )ds + aE(M t ) + aW t - M t , t> 0, (1.2) 
Jo 

where Xq < 1 almost surely, a G M, cr > 0, Wt is a standard Brownian motion in 
M and b : R — > R is Lipschitz continuous. The jumps, or resets, of the system are 
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described by 



k > 1 



where (rk)k ^ 1 stands for the sequence of hitting times of 1 by the process {X t )t ^ o- 
That is (M t )t > o counts the number of times X t hits the threshold before time t, so 
that E(Mt) denotes the expected number of times the threshold is reached before t. 
This equation corresponds to what we would expect when the limit N — > oo is taken 



in (1.1 ) when the above simplifications are made, where the resets are now included 
in the equation by the presence of the counting process (M t ) t . It is indeed the 
case that the associated finite system of N neurons converges to this nonlinear limit 
equation (again this is shown in our forthcoming manuscript [5J). 



In fact a less general version of equation (1.2) has been rigorously studied from 
the PDE viewpoint before, though not completely satisfactorily. If one applies 
the generalized ltd formula for jump processes to f(X t ) for a suitable function /, 
then in the standard way one arrives at the Fokker-Planck equation for the density 
p(t,y)dy = F(X t edy): 

d t p(t,y) + d y [(%) + ae'(t))p{t,y)) - -d 2 yy p{t,x) = 5 Q (y - V R )e'(t), y < 1, 

where e(t) = E(M t ), subject to p(t, 1) = 0, p(t, -oo) = 0, p(0,y)dy = P(X G dy) 
and where So is the Dirac measure at 0. Moreover, the condition that p(t, y) must 
remain a probability density translates into the fact that we must have 

e'(t) = j^{M t ) = -±d y p{l,t), Vt > 0. 

Now, in the case when b(x) = —x, this nonlinear Fokker-Planck equation is exactly 
the one studied in [3] and [4]. Therein, the authors define a solution to be one in 
which e'(t) remains bounded, since the explosion of e'(t) corresponds to the situation 
where a large proportion of the neurons in the network all spike at the same time. 
Given this definition, the equation is shown to be uniquely solvable in the cases 
a = and a < (see Theorem 1.1 in [4|), the latter one being referred to as 
'self-inhibitory' in neuroscience. In the so-called 'self-excitatory' framework, i.e. 
for a > 0, existence of a solution for all time has been left open. Instead, very 
interestingly, a negative result was established in j3], stating that, for any a > 0, it 
is possible to find an initial probability distribution P(Xo G dy) such that there is a 
blow-up, i.e. such that e'(t) = oo in finite time. As a remarkable point, it says that 
solvability in the long run may fail for possibly small values of a. On the contrary, 
it does not come as a surprise that solutions may explode in a finite time for the 



prescribed values of Vr and Vp when a ^ 1: indeed, if Np particles in ( 1.1 ) reach the 
threshold Vp = 1 between some times ti and t% for some proportion p G (0, 1], then 
all the particles perform a jump of size greater than ap(N — 1)/N within the same 
period, which reads as a size greater than ap in the ideal case N = oo corresponding 



to (1.2); in particular, when ap ^ 1 and the reset is fixed as 0, particles keep on 



crossing the threshold an infinite number of times between t\ and ti- Such a picture 
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indicates that the system enters a sort of 'super-self-excitatory' regime when a is 
above Vp — Vr, which matches 1 with our convention. 

Here we thus investigate the case a G (0, 1). We then contribute the associated 
positive and until now unknown result, by proving that, given a starting point 
X Q = xq, we can find an explicit a small enough so that there is an unique solution 
to (1.2) (and hence to the associated Fokker-Planck equation) for all time, such that 
e'(t) = [d / dt]E(M t ) remains bounded on any compact subset (see Theorem 5.2). 
The surprising difficulty of this result is reflected in the rather involved nature of 
our proofs. 

Our result is for a general Lipschitz function b, but there are two important spe- 
cific cases that we study further: the Brownian case when 6 = and the Ornstein- 
Uhlenbeck case when b(x) = —Xx, A ^ 0. The Ornstein-Uhlenbeck case is most 
relevant to neuroscience since it is often used in the literature, but surprising diffi- 
culties remain in the purely Brownian case. In both these cases we are able to give 
an explicit a depending on the starting point Xq such that (1.2) has a solution for 
all a < Oq. However, it is interesting to note that our explicit values do not appear 
to be optimal, in the sense that from our numerical simulations there appears to 
be a real 'gap' between our positive result and the negative result of |3]. Precisely, 
our simulations show that for a given xq there exist classical solutions for all time 
for a bigger than our explicit ato, while there exist solutions that blow-up that do 
not satisfy the conditions of (3]. Thus a very interesting, and perhaps very difficult 
question is to determine for a given initial starting point the true critical value a c 
such that for a < a c (1.2) has an unique global solution. We discuss this further at 
the end of the paper. 

Equation (1.2) can be thought as of McKean-Vlasov-type, since the process X t 
depends on the distribution of the solution itself. However, it is highly non-standard, 
since it actually depends on the distribution of the hitting times of the threshold by 
the solution. This renders the traditional approaches to McKean-Vlasov equations 
and propagation of chaos, such as those presented by Sznitman in [20], inapplicable, 
because we have no a priori smoothness on the law of the hitting times. Thus our 
results are also new in this context. 

The general structure of the proof is quite typical of the methods used to investi- 
gate well-posedness of nonlinear systems, and especially well-posedness of Markovian 
stochastic differential equations (or equivalently of partial differential equations of 
second order) involving some non-trivial nonlinearity. Precisely, the first point is to 
tackle unique solvability in small time: when the parameter a is (strictly) less than 
1 and the mass of the initial condition decays linearly at the threshold, it is proven 
that the system induces a natural contraction in a well-chosen space provided the 
time duration is small enough. In this framework, the specific notion of solution at 
hand plays a crucial role as it defines the right space for the contraction. Below, 
solutions are sought in such a way that the mapping t i— > E(M t ) is continuously 
different iable. The proof is at the intersection between probability and PDEs. The 
second stage is then to extend existence and uniqueness from short to long times. 
Exactly as for PDEs of quasilinear type, the point is to prove that some key quantity 
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is preserved as time goes by. Here, we prove that the system cannot accumulate too 
much mass in the right vicinity of 1. Equivalently, this amounts to showing that 
the Lipschitz constant of the mapping t h-> K(M t ) cannot blow-up in a finite time. 
This is where the condition a small enough comes in: when a is small enough, we 
manage to give some estimates for the density of X t in the neighbourhood of 1, 
the critical value of a explicitly depending upon the available bound of the density. 
Generally speaking, we make use of standard Gaussian estimates of Aronson type of 
the density; in the cases b = and b(x) = —Xx, for some A > 0, the bounds 

are explicit enough to derive a numerical ao up to which unique solvability holds. 
Unfortunately, the estimates we use are rather poor as they mostly neglect the right 
behavior of the density of X t at the boundary, thus yielding a non-optimal value. 
Anyhow, they serve as a starting point for proving a refined estimate of the gradient 
of the density at the boundary: this is the required ingredient for proving that, at 
any time t, the mass of X t decays linearly in the neighbourhood of 1, uniformly 
in t in compact sets, and thus to apply the existence and uniqueness argument in 
small time. Such a refined estimate is proven on any interval where a solution does 
exist: below, it is thus referred to as an a priori estimate. Coupled with the result 
of existence and uniqueness in small time, it permits to prove by induction that 
existence and uniqueness hold on any finite interval and thus on the whole of [0, oo). 

The layout of the paper is as follows. In the first section below we make the 
necessary definitions and describe the setting in which we will be working. Section 



3jis devoted to proving the existence and uniqueness of a solution to (1.2) in small 
time. Long-time a priori estimates are proven in Section |4j This section contains 
the crux of our argument, since we here show that, for a small enough, if there exists 
a solution, it does not exhibit the blow-up property highlighted in [3]. In Section [5] 
the final step is made, whereby we show that under the a priori Lipschitz bound, an 



unique solution exists to (1.2) that does not blow-up at any time. We conclude the 
paper in Section [6] with a discussion of the critical values of a in both the Brownian 
and the Ornstein-Uhlenbeck cases. 

2. Setting 

As stated in the introduction, we are interested in solutions to the nonlinear 
McKean-Vlasov-type SDE 

X t = X + [ b{X s )ds + aE{M t ) + W t -M t , t^O, (2.1) 



o 



where X < 1 almost surely, a £ (0, 1) and W t is a standard Brownian motion with 
respect to a filtration (J-"t)t > o satisfying the usual conditions, with values in R and 
starting from 0. The jumps, or resets, of the system are described by 

M t =^l M (T k ), (2.2) 

k > 1 

where the sequence of stopping times (rfc)* > o are defined by tq = and 

T k = M{t > r fc _x : AV ^ 1}, k > 1. (2.3) 
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We assume that b : (—00, 1] — > R is Lipschitz continuous such that 

\b(x)\ < A(|x| + 1), \b(x) -b(y)\ ^ K\x - y\, Vx, y e (-00, 1]. 



Remark 2.1. Although we take the intensity of the noise in (2.1) to be 1 rather than 
a general a > as stated in the introduction, this results in no loss of generality. 
Indeed, by setting u = t/o~ 2 we can rewrite (2.1) as 

X a2u = X + b(X s )ds + oE(M$ u ) + W a2u - M* 2u , u^O, 
Jo 

for an arbitrary o > 0, where we have emphasized the dependence of the process M 
on X by writing M = M x . Now, by Brownian scaling we have that W a 2 u = aW u 
where W is another Brownian motion starting from 0. Set X u = X a 2 u . Then, 
since the number of times the process X hits the threshold before time a 2 u is equal 
to the number of time the process X hits the threshold before time u, we have that 
M£ u = M*. Hence it follows that 

fG^U _ _ 

X U = X + / b(X s )ds + aE(M* ) + oW u - M* , u^O. 
Jo 

By a change of variable we then have 

X U = X + / a 2 b(X s )ds + aE(M u x ) + aW u - M* , u ^ 0. 



Thus, if we modify the function b accordingly, X satisfies exactly the equation of 
interest, but now with an arbitrary noise intensity. 

As discussed in the introduction, the key point is to look for a solution for which 
t 1 — y K(M t ) is continuously differentiable, which would correspond to a solution 
whose derivative does not exhibit a finite time blow-up. To this end we introduce 
the following objects: 

• The space C: For T,A^0, let 

C(T,A) = \e eC^O,! 7 ] : e(0) = 0, e(s) < e{t) Vs ^ t, sup |e'(t)|<Al 

I s£ t sg T ) 

where as usual C [0,T"] denotes the space of continuously differentiable functions on 
[0,T]. 

• The map Y: For e € £(T, A) define 

r(e)(t):=E(M t e ) (2.4) 

where 

Mt= J> M (r fc e ) (2.5) 

k > 1 

and = inf{t > r^_ 1 : X^_ ^ 1}, r = and (X t e , M t e ) is a solution to 



X e t = X + f b(X e s )ds + ae(t) + W t - M 4 e , t ^ 0, X < 1. 
Jo 



(2.6) 
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Note that the solution is well defined as b is Lipschitz continuous. The map T is 
thus defined as a map from £(T, A) into the set of non-decreasing functions on [0, T\. 
It thus depends on A as its domain of definition depends on A; for this reason, it 
should be denoted by Y A . Anyhow, since the family {T a )a > o is consistent in the 
sense that, for any A' ^ A, the restriction of T A to C(T, A') coincides with T A , we 
can get rid of the superscript A in T A and then use the simpler notation T. 

The idea is to look for a fixed point of T in a certain space, which means that we 
are to find a space that is kept invariant by T. With this in mind, the following a 
priori stability result for F will be useful, the proof of which we leave until the end 
of the section. 

Proposition 2.2. Given T > 0, e G C(T,A) and some real a ^ E[(Xo)+] (where 
(x) + indicates the positive part of x G M.), it holds that 



with 



(Vt G [0,T], e(t) < g(t)) =}> (Vt G [0,T], T(e)(t) < a(t)). 
a + (4 + AT 1/2 )t 1/2 



(2.7) 



cxp 



2At 



1 — a \1 — a 

With this in mind, we thus also introduce the space 7i, followed by our definition 



of a solution to (2.1): 



The space H: For an initial value Xq < 1, let 

H(T,A,X ) = {e G C(T,A) : e(t) < C T {Vt + E|X |)e e '} 



where CV and # are chosen according to Proposition 2.2, that is 

4 + AT 1 / 2 „ 2A 



T ■- 



1 — a 



9 :-. 



1 — a 



Equip "H(T, A,X ) with the norm ||e|| w(T)Ai x ) = ||e|| 00) T+||e'|| D0) r. Then "H(T, A, X ) 
is a complete metric space, since it is a closed subspace of C 1 [0, T\. We then define: 



Definition 2.3 (Solution to (2.1 )). The coupled process (X t , M t )o <g t < r wi/Z sazd 
6e a solution to (2.1) nj> ■untzZ timeT if (M t )o^t^T satisfies (2.2) subject to (2.3), 



([0,T] 9 t i-> E(Mtj) G H(A,r,I ) /or some A > 0, and (X t ) 0<t <T *s « s£ron# 
solution of ( 2.1[ ) up nnizZ izme T. 



With the above definitions, it is clear that, for any A ^ 0, any fixed point of T 
that belongs to H(A,T, Xq) provides a solution to ( |2.1 ), and conversely that, for 
any solution (X t , M 4 ) < t < t of (2.1 ), ([0, T] 9 t 1— >■ E(M 4 )) must be a fixed point of 
T in TL(A, T, X ) for A as in Definition 2.3 



2.1. Proof of Proposition 2.2 



Proof of Proposition 2.2: Fix T > 0. We first note that we may write 



M t e = supL«)J, Zt = Xt + M^ tG[0,T], 



(2- 



s < t 
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Indeed, one can see that for t G [t^,t^ +1 ), k ^ 0, 



sup [(Z e s ) + \ = max ( sup [(X e s + j) \ , sup [(X e s + k) + \ ) 

S<t 0<J<fc-l ^ a6 [ T /,T/ +1 ) S<=[T*,t) J 



max (j + 1, AO = M, 

< j < fc_i v 7 



t > 



using the fact that X 4 e < 1 for all t ^ 0. 

Then, given t G [0,T] such that Z e t ^ 0, let p e : = sup{s G [0,t] : Z% < 0} 
(sup0 = 0). Pay attention that p e is not a stopping time. Then, for s G [p e , £], 

\b(x:)\ < a(i + |x:d < a(i + + m;) = a(i + (z e s ) + + mi). (2.9) 



By (2.8), we know that M e s ^ sup < r < 8 {Zf) + . Therefore, 



\b{X e s )\ < A ( 1 + 2 sup (Z r e 

< r < s 



By (2.6), we obtain: 



Z*^Z e pe +A 1 + 2 sup (Z r e )+ rfs + ae(t) +W t -W. 



sir 



(2.10) 



If p e > 0, then Zp B = as, obviously, (Z^) ^ s ^ y is a continuous process. If p e = 0, 
then Zpe = Xq ^ 0. Therefore, 

(Z t e )+ < (X )+ + A / fl + 2 sup (Z T e )+)rfs + ae(t) + 2 sup (2.11) 

JO \ OsCrsJs / s: s s£ i 

Obviously, the above inequality still holds if Zf ^ 0. Therefore, taking the supre- 
mum in the left-hand side, taking the expectation and applying Doob's maximal 
inequality for martingales, 



E 



sup {Zt) + 

< s < t 









sup (Z e r ) + 




.0 < r < s 



which shows that E[sup < t < T (Zf) + ] is finite. In particular, if 



ds + ae(t)+At 1/2 , 
(2.12) 



VtG[0,T], e(t)^g(t), with (/(*):= 



E[(X )+] + (4 + AT 1 / 2 )* 1 / 2 / 2At 



1 — a 



cxp 



1 — a 



and i? e stands for the deterministic hitting time: 



i? e := inf h G [0,T] : E 



sup (Z c s 

€ s < t 



>g(t)\ (inf0 = +oo) 



AN INTEGRATE-AND-FIRE MODEL OF MCKEAN- VLASOV TYPE 



9 



then, for any t G [0, R e AT], 



E 



sup {Zl 

€ s <. t 



< E[(X ). 



< (E[(X ) + 
= (1 - a)g(t) + ag(t) 



A / (l + 2g(s))ds + ag(t) + 4t 1/2 
•/ o 

- (4 + AT 1 ' 2 )?' 2 ) 



2A 



2As 



a: 



ds 



+ as(*) 



t s; Tl^t J+J> 



Now, by the continuity of the paths of Z e and by the finiteness of E [sup 
we deduce that E[sup ^ s ^ t (Z e s ) + ] is continuous in t. Therefore, if R e ^ T, then 
E[sup 0< s< Re(Z%) + ] must be equal to g(R e ), but, by the above inequalities, this 
sounds as a contradiction. By (2.8), this proves the announced bound when a in the 

] . The proof is completed in the same way when 

□ 



statement exactly matches E[(X ) 
a>E[{X ) + ). 

Remark 2.4. In the case when Xq = 0, we have that T exhibits the same stability 
on functions e such that e(t) ^ Cry/ie 91 . 

3. Existence and uniqueness in small time 

The main result of this section is the following: 

Theorem 3.1. Suppose there exist (3,e > such that F(Xq G dx) ^ /3(1 — x)dx 
for any x G (1 — e, 1] and that the resulting density of the law of Xq on the interval 
(1 — e, 1] is differentiate at point 1. Then there exist constants A\ ^ 1 and T\ G 
(0,1], depending upon X through (5, e and E|X | (which is assumed to be finite) 
only, such that 

T(H(A u T l7 X )) CH(A 1 ,T U X ). 
Moreover, for all ei,e 2 G %{A 1 , T l5 X ), 

1 



|r(e 1 )-r(e 2 )|| 



H(A lt T lt Xo) 



ei - e 2 



H(Ai,Ti,X ) 



Hence there exists an unique fixed point of the restriction of V to H(Ai,Ti,Xq) 
which provides a solution to ( 2.1[ ) according to Definition 2.3 up until time T\. 



Remark 3.2. The assumption F(Xq G dx) ^ (3(l — x)dx, forx G (1 — e, 1], says that 
the distribution of X has a density on the interval (1 — e, 1], and that the density is 
bounded by (3(1 — x). This is a key assumption to prove Thereom 3.1 since it prevents 



the hitting time of 1 by X from having a Dirac mass at the initial time 0. In the 
proof, it permits to show that T(e) is Lipschitz continuous when e itself is Lipschitz 
continuous. By contrast, we feel that the proof could be written without requiring 
the density of Xq at point 1 to be differentiate. Here the differentiability of the 
density forces T(e) to be continuously differentiable at time when e is continuously 
differentiable, as specified in the definition of C(T,A). This proves to be really 
convenient in the whole analysis, but does not affect the final result in Section [5[ 
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3.1. PDEs. Given e G C(T,A), we have 

r(e)(t) = E(M t e ) = G M W = *) P W G rfs) + P(r 1 e < t), 

where P(t^ G ds) is a convenient abuse of notation for denoting the law of rf. and 
(S(R) 9 A i— ?• P(r|, x G A|r| = s) stands for the conditional law of r^ +1 given r| = s. 
Here 13 (M.) is the Borel cr-algebra on M. Observing that the solution X e to (2.6) is 
a Markov process, we may write 

r(e)(t) = E(M t e ) = V f ¥(rf s <: t - s\xf = 0)P(r A e G ds) + P(r 1 e < t), (3.1) 

where e" 3 stands for the mapping ([0, T — s) 3 t i— > e(t + s) — e(s)) G £(T — s,A). 
In order to understand the behaviour of T, it will become clear that one must study 
the marginal density of the killed process 

Note here we have used the notation of the Radon-Nikodym derivative to denote 
the density of the killed process with respect to the Lebesgue measure, which makes 
sense thanks to the following Lemma: 

Lemma 3.3. Given an initial condition Xq satisfying the same assumptions as 
in the statement of Theorem 3.1 together with e G £(T,A), consider the solution 
(Xt)o<t<T of the SDE 

d X t = b{ X t)dt + ae\t)dt + dW t , t G [0, T] ; X o = X , (3.2) 

together with the stopping time 

n = w£{te[0,T\: X t>l}- 
Then, for any t G (0, T], the measure 

B((-oo, 1]) 3 A ^ F( Xt G A, t < n) 

has a density with respect to the Lebesgue measure. Here, B((— oo, 1]) denotes the 
Borel subsets of (— oo, 1]. Pay attention that the density is not normalized. Denoting 
the density by 

p{t,y) = ^F( Xt Edy,t< n), t G [0,T], y < 1, (3.3) 

it is continuous and continuously differentiable in y on (0, T] x (— oo, 1] and admits 
Sobolev derivatives of order 1 in t and of order 2 in y in any L q , q ^ 1, on any 



compact subset of (0, T] x (— oo, 1). When f3 — in the statement of Theorem 3.1 



i.e. Xq ^ 1 — e a.s., it is actually continuous and continuously differentiable in y on 
any compact subset o/([0,T] x (— oo, 1]) \ ({0} x (— oo, 1 — e]). Moreover, p satisfies 
the Kolmogorov equation: 

d t p(t,y)+d y [(b(y)+ae'(t))p(t,y)] - l -d 2 yy p{t,y) = 0, t G (0,T], y < 1, (3.4) 
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with the Dirichlet boundary condition p(t, 1) = and the measure-valued initial 
condition p(0,y)dy = P(X G dy), both p(t,y) and d y p(t,y) decaying to as y — >• 
— oo. Finally, T\ has a density on [0, T], given by 

d 



dt 



(ri < t) 



-d yP (t,l), te[0,T], 



(3.5) 



the mapping [0, T] 3 t i— > d y p(t, 1) being continuous and its supremum norm being 
bounded in terms ofT, a, A, (5 and b only. 



Lemma |3.3| is quite standard. Anyhow, it is difficult to find a complete proof of 
it under our assumptions. For this reason, we provide a complete proof in appendix 
(Section [7}. We can apply the above result to (13 . 1 h to obtain: 

Proposition 3.4. Given an initial condition Xq satisfying the same assumption as 
in the statement of Theorem 3.1 together with e G C(T,A), the mapping [0, T] 3 
1 1 — y r(e)(i) is continuously differentiable. Moreover, 



d 
dt 



[r(e)] (t) = - jf \d y pf< s \t- Sl 1)A [r(e)] (a)ds- i^p.(t, 1), t G [0, T], (3.6) 



where p e represents the density of the process X e killed at 1 and p^e'^ represents the 
density of the process X e9s killed at 1 with Xf = 0. 

Proof. We first check that T(e) is Lipschitz continuous on [0, T]. Considering a finite 
difference in (3.1 ) and using ( |3.5[ ), we get, for t, t + h G [0, T], 

r(e)(t + h)- r(e)(t) = ^ [ t+H F(rf s ^t + h- s\xf = 0)P(r fc e G ds) 



k > 1 
■t+h 



t pt+h-s 



d y p^ s \r,l)drF(r^ G tfe) 



(3.7) 



Jt-s 



By Lemma 3J3 we can handle the two last terms in the above to find a constant 
C > such that 

r( e )(t + /i) - r( e )(t) < Yl / p ( r i* s ^ * + h - s \ x o a = °) F ( T k G ds ) 

+ C7i(l + r(e)(T)), 



the last term in the right-hand side being finite thanks to Proposition 2J2 Moreover, 
from the continuous differentiability of e, we deduce that 



lim sup P(r 1 e ^ h\Xf = 0) 

h \°0 € s <T-h 



lim sup P 

h \°0 <S s < T-h 



sup ^ l|X e " 

< r < h 



(3.8) 
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where Z e s is given by (2.8), the second line following from (2.11). Therefore, there 
exists a mapping rj : M.+ — > R + matching at and continuous at such that 

T(e)(t + h)-T(e)(t) ^ rj(h) [V(e)(t + h) - V{e){t)} + Ch(l + T(e)(T)) . 

Choosing h small enough, Lipschitz continuity easily follows. 

As a consequence, we can divide both sides of (3.7) by h and then let h tend to 
0. By (3.8), we have for a given t G [0,T), 

-t+h 



lim hT 1 V / F(rf s ^t + h-s\Xf = 0)P(r fc e G ds) 
hS, ° k>i Jt 

r(e)(f + /i)-r( C )(f) 



^ lim 

7i\Q 



sup P(rf ^ /i|X e ° a =0)- 

< s < T-h 



h 



0. 



Handling the second term in (3.7) by Lemma 3.3 and using the Lebesgue Dominated 
Convergence Theorem, we deduce that 



d 
dt 



F(e)(t) 



By Lemma 3.3, we know that dyP^' s \-, 1) and d y p e (-, 1) are continuous (in t). This 
proves that (d/dt)T(e) is continuous as well. 



Formula (3.6) follows from the relationship 

-t 



(3.9) 



□ 



3.2. Density estimates in short time. Suppose that X = x < 1. Following 
the parametrix method in P, Chapter 1], one can show that the solution to (3.4) 
(then denoted by p^° ) is given by 



p?(t,y) = q{t,x ,y) 



(ae'(s) + 6(1)) d z pf (s, z)q(t - s, z, y)dzds } (3.10) 



for t G (0, T] and y < 1, where q(t,x,y) is a solution to the PDE 

d t q(t, x, y) = (l/2)d 2 yy q(t, x, y) - d y [(%) - b(l))q(t, x, y)] 
q(t, x,l) =0 
Q@,x,y) =5 (x-y) 



(3.11) 



on [0, T] x (— oo, 1] x (— oo, 1]. We will need the following bound on the behaviour 
of q. 



Proposition 3.5. For any T > there exists a constant B? > such that 



B 



\d y q(t,x,y)\ < -^exp 



Iff -2/1 
B T t 



(3.12) 
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and 

\d y q{t,x,y)\ ^ ^exp (- |g ' ~ B ff ^j , (3.13) 
for all t G [0,T] ; x, y ^ 1, where 

Q = x + [\g)dr, tGM, (3.14) 
Jo 

and b(v) = b(v) — 6(1) if v < 1 while b(v) = —6(2 — v) if v ^ 1. 

Proof. Let (X t ) t ^ o denote the solution of the SDE: 

dX t = b(X t )dt + dW t . (3.15) 

We can then see that q(t,x,y) is the transition density of X t killed at 1 (equation 
(3.11 ) is the Fokker-Planck equation for the killed process, as shown in Lemma 3.3). 

Note that 6(1 + v) — —6(1 — v) for v G K, so that 6 is odd with respect to 1. 
Therefore, when initialised at Xq = 1 the process (X t — l)o ^t^T satisfies 

X t - 1 = / b(X s - 1 + l)ds + W U t G [0, T] 
Jo 

and the process (1 — X t )o ^t^T satisfies 

1 - X t = - I b(X s )ds -W t = [ 6(2 - X s )ds - W t 
Jo Jo 

= [ b(l-X s + l)ds-W t , 
Jo 



which is the same equation. Hence (1 — X t ) ^t^r and (Xt — l)o < t < r have the 
same distribution when X = 1. Therefore the reflection principle (see [T?j Section 
1.13]) applies and 

q(t,x,y) = q(t,x,y) -q(t,x,2-y), te[0,T], x,y^l, (3.16) 

where q(t, x, y) is the transition density of the non-killed process. 

We thus need to estimate the transition density q of the non-killed process. By 
Theorem 1.1 of [6] it follows that 

g(M,j/)<%exp("- |y ~f |2 y t>0, x,yeR, (3.17) 



y/t V B Tt 
where 

e? = x + f b(&ds, t > o. 

JO 

We claim that a similar bound holds for \d y q(t,x,y)\, i.e. 



\d y q(t,x,y)\ < ^exp (~ V ^ ) , t > 0, x,t,6i (3.18) 
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Before we prove (3.18), we show that it implies (3.12) and (3.13). By (3.16), we 
deduce that 



\d y q(t,x,y)\ ^ 



B 



r 



cxp 



|y-ff 12 

B T t 



+ exp 



y - it 



x\2 



B T t 



t > 0, x,y ^ 1. 



Since 6(1) = 0, we have = 1 for any t ^ 0. Therefore, by the comparison principle 
for ODEs, it must hold £f ^ = 1 for any t ^ and x ^ 1. As a consequence, for 
t > and x,y ^ 1, 

|y-e?l = l(y-i)-(€f-i)l<|y-i| + l€f-i| 

= 1-^ + 1-^ = 2-^^ = 12-2/-^. 

Therefore, the largest exponential term in the above bound for \d y q(t,x,y)\ is the 
first one. Equation (3.12) follows directly (adjusting Bt as necessary). Now, by 
Gronwall's Lemma and thanks to the Lipschitz property of b, we have: 



eM-Ks)\g - y| < \Cts - e s \ < exp(Ksm ~ v\ 



(3.19) 



for all s e [0,t], so that, choosing s — t, we get (3.13). 

We now turn to the proof of (3.18). Without any loss of generality, we can 
assume b to be a twice continuously differentiable function with a bounded second- 
order derivative. Indeed, if we can prove that, in such a case, (3.18) holds with 
respect to a constant Bt that does not refer to the second-order differentiability 
of b, then (3.18) holds in the original setting as well by a standard mollification 
argument. Since q is the transition probability of the solution to (3.15), it satisfies 
the Fokker-Planck equation 

d t q(t,x,y) = -dy y q(t,x,y) -b(y)d y q(t,x,y) -b'(y)q(t,x,y), 



for t G [0, T] and x, y G R. Define (Xt)t > o to be the solution of the SDE 

-b(Y t y )dt + dW t , 



dY t y 

1 



y- 



(3.20) 



By Theorem 39 [TBI p. 312], we know that, a.s., the mapping [0, +oo) X R 3 (t, y) i— > 
Yf is continuously differentiable with respect to y and that 



d y Y£ = exp 



{-fm 



)du . 



(3.21) 



Moreover, by Lemma 3.3, we know that q is regular enough in order to apply the 
Ito-Krylov formula (see [T2l Section 11.10]) to [0,t/2] 3 s \— > q — s,x,Y^j, which 
yields 

-t/2 



q(t,x,y) = E q it/2,x,Y t y /2 



E 



V(Y?)q[t-8,x,Y? 



ds. 
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By the Malliavin-Bismut-Elworthy formula (see, for example, Theorem 2.1 of [7]) 



d y q(t,x,y) = -E 



q(t/2,x,Y? /2 ) 



t/2 



d y Y s y dW s 



t/2 



-E 



bf(yv)q(t-s,x,Yy) I d y Y r ydw r 



ds. (3.22) 



By (3.21), ISyY^I ^ exp(Ks). Thus using (3.17) we can compute for any < s ^ t/2 
(where the constant Bt changes from line to line below) 



I(s,t) :-- 



S L 

Bt 



-E q(t-s,x,Y s y ) I d y Y r y dW r 
E[f(t-s,x,Y s y ) 

1 



" 1/2 






Hi 







S L 

Bt 



1/2 



Bt 



t - s 

1 

t - s 



E exp 



exp 



yy _ cx 12 

1 s St-sl 



1/2 



1/2 



5 r (t - s) 



B T (t-s) 



P(Y/ e dz) 



1/2 



By an estimate similar to (|3.17|) but for the density of Y, we deduce 
Bt 



I(s,t) ^ 



s [(t - s)^/s 



exp 



Bt 



exp 



l~ t& 1 2 

B T (t - s) 

\tx _ tV 12 

B~r~t 



exp 

1/2 



dz 



1/2 



where, for the first line, we have used the fact the flow associated with the ODE 
driven by —b is nothing but the backward flow associated with the ODE driven by 
b, i.e. (£* t )t ;s o satisfies 



£ t = x- b(Cl s )ds, t > 0, 
Jo 

and, to pass from the first to the second line, we have recognised a Gaussian convo- 
lution. By (3.19), we deduce that 
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from which ( |3.22[ ) yields 

\d y q{t,x,y)\ < — exp 



+ B T exp 



B T t 

r t -y\ 2 \ f t/2 i 



B T t 



-,1/2 



y/t^sy/t. 



ds 



Bx 
^— exp 



( \it-y 


to 


+ 7f exp 


( ier-y 


to 


V B T t ) 





(is. 



This proves (3.18). 



□ 



Lemma 3.6. Let T > 0. Suppose X = xq < 1 and e G £(T, A). T6en there exists 
a constant K\(T) (independent of A and Xq) which increases with T such that 

sup [vG\\d y p*°(s,.)-d y q{s,x ,.)\\J <(aA + |6(l)|)«i(r), 

< s < t 

/or a// non-negative t ^ min{[(o;yl+ |6(l)|)/ti(T)] -2 , T}. 
Proof. For notational sake, define 

i^°(t,2/) :=d„^ (t,y)-3 y g(t,x ,2/). 



By (3.10), we see that 
\F*°(t,y)\ ^(«A +|6(1)|) 

■t rl 



t pi 



\d z q(s,x ,z)d y q(t - s,z,y)\dzds 



J-oo 



(3.23) 



+ (aA+\b(l)\) f I \d z p x e °(s } z) - d z q(s,x ,z)\\d y q(t - s,z,y)\dzds. 

JO J-oo 

Thus defining A := aA + |6(1)| and allowing By to increase as necessary from line 



to line below, by Proposition |3.5| we see that for t ^ T 

\F?{t,y)\ 



t pi 



< AB T / / s" 1 (t-s)~ 1 exp 



J-oo 



( 






exp 


( 




\ B T s J 


\ B T (t - s) 



dzds 



+ AB T / (t-s)- 1 \\F*°{s 



1 I II oo 



exp 



(t- s )l 



dzds. 



B T (t - s) 

We can recognise the first term in the above as a Gaussian convolution. Thus 



i*r(*,v)i< 



A3 



T 



\/t Jo y/sy/t - S 



exp 



AB T 
1 



F5°(s, OIL ds 



< AB T ( — = exp 



ie t -*oi 2 



(3.24) 



B?t 



+ sup [v5P7°0v)IU 

< s < t 
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where again we have used (3.19). Therefore, by multiplying both sides of the in- 
equality by y/i, we deduce 

sup [Villi^^OlU^ABr + V^ABr sup [sfs \\F*°{s, .)||J , 



^ s < t 



< s < t 



We can conclude that, for t G [0,T] such that i < [2AB T ]~ 2 = [2{aA+\b(l)\)B T ] 



we have that 



sup [ v ^\\d y pt°(s,.)-d y q(s,x (h .)\\ oo \ ^2AB T . 



< s < t 



(3.25) 
□ 



Proposition 3.7. Lei T > 0. Suppose Xq — xq < 1, and that e G C(T,A). Then 
there exists a constant K^CO (independent of A and Xq) which increases with T such 
that 

1 ( 1^0 _ y\ 2 * 

\d y p:°(t,y)-d y q(t,x ,y)\ s= K 2 (T)(aA + \b(l)\)-^ex V I - ^ 

/or a// £ ^ min{[(ayl + |5(l)|)/t2(T)]~ 2 , T} and y ^ 1, where (Ct°)t^o ^ s as i> n 
Proposition \3.5\ 



Proof. Throughout the proof, we will use the fact the mapping 

exp(x) — 1 



if : K 3 x h-> 



,r 



(V(0) = 1), 



(3.26) 



is non-decreasing. We will also use the same notation as in the proof of Lemma 3.6 

F?(t,y) :=d y p: o (t,y)-d y q(t,x ,y). 



First Step. As noted in (3.23), we have for t G [0,T] 



\F*°(t,y)\^A [ [ \d z q(s,x ,z)d y q(t-s,z,y)\dzds 

JO J -co 

+ l[ I \F:°(t 1 ,y 1 )\\d y q(t-t 1 ,y 1 ,y)\dy 1 dt 1 

JO J-oo 



(3.27) 



where we have denoted A = olA + |6(1)|. The fact that if given by (3.26) is non- 
decreasing says that, for t G [0, T], 

1 - exp(-2ii~T) < 1 - exp(-2Kt) < ^ < exp(2Kt) - 1 < exp(2fs:T) - 1 



2KT 



2Kt 



2Kt 



2KT 



from which, together with Proposition 3.5, we deduce there exists > such that 

tv 1 2 \ 

(3.28) 



Mt,x,y)\ < IZ ^exp (- *^* 



\d y q(t,x,y)\ ^ 



Bn 



2Kt 



exp 



2B T (e 2Kt - 1) 



(3.29) 
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Unlike in the previous lemma, Bt is now fixed. We then introduce the kernel 



1 



K(t,x,y) = < 



1 



{ y/B T [l-e 



2Kt^ 



exp 
exp 



\it-y\ 



2B T (e 2Kt - 1) 
2B T (1- e 2Kt ) 



if t > 0, 
if t < 0. 



(3.30) 



For < s < t ^ T, /C satisfies the Gaussian convolution property: 

JC(s, x, •) <g) ^C(— (t — s), -, i/) := ^ /C(s, x, z)/C(— (i — s), z,y)dz 

: exp 1 



y/B T (e 2Ks - e - 2K ^- s )) \ 2B T (e 2Ks - e ^ 2K ^- s )) 
By fl3~T9| ), we know that |f£ - y\ < e ^ (t - s) |^ - £^ (t _J for all < s < t. Therefore, 



/C(s, re, -)® £(-(*- s),.,y) < 



exp 



«C V2^e KT IC(t,x,y) 



l€f-y| s 



2Br( e 2«-l); ( 3.3i) 



By ( pQ7| ), ( |08| ) and ( pQ9[ ), we deduce that, for i 6 [0,T], 



|F^(t,|/)|^^e KJ /C(t,Xo,i/) 



o Ve 2As - Wl - e- 2K ^~ s ) 

' ,l \F e X0 (ti,yi)\ 



ds 




(3.32) 



>o j-oo Vl - e- 2i ^-'i) 
Using (3.26), we notice that 

V2KT 



)C(-(t - h),y l ,y)dyxdt 1 



1 



Ve 2Ks - - e~ 2K ^- 



,ds 



1 - exp(-2KT) Jo V2Ks^2K(t - s 



-.ds 



T 



y/2K(l - exp(-2KT)) Jo y/sy/t - s 

-. Or, 



ds 



y/2K(l - exp(-2KT)) 
using § Q \\fs\/i — s)~ 1 ds = L(\/u\/l — u)~ l du = n ^ 4. Finally, (3.32) yields 
|F^(t,y)|^icVi^/C(t,xo,y) 



(3.33) 



iq J-oo Vl - e- 2 ^(*-*i) 
with <5 T := v / 2^C T exp( J fs:T). 



' '"' IF" ^ will ( 3 - 34 ) 
1 ' ' ! -" /J " K(-(t-t 1 ),y h y)dy 1 dt h 
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Second Step. We now prove by induction, that for any N ^ 0, 

TV 

\F?(t,y)\ ^ AC T B^(t,x ,y)Y,(^^C T B 3 T /2 Vt) 1 + R N+1 {t,y), (3.35) 



i=0 



where 

R N (t,y) = {AB 



3/2x/V 
T ) 



N-l 



dt 



N 



s: t N < ... < t 



i=0 



-2if(ti-* i+ i) 



„ N-l 

/ dy N ...dy^F^^N.y^] TT /C(-(*j - t i+1 ),y i+1 ,yi), 

J(-oo,l]» to 

with the convention t = t and y = 2/- 111 the first step, we established (3.35) when 
N = 0. Assume now that, for some N ^ 0, (3.35) holds for any t G [0, T] and y ^ 1. 
Then, plugging the induction assumption at rank iV into (3.34), we get: 

\F:°(t,y)\ 

^AC T B*JC(t,x ,y) 



/ N i\ r 

+ AC T B Z T (VKAC T BrVty) AB^ 2 J 



\t r3 /2 [' /C(ti,g ,-)(8)/C(-(t-ti),-,y) 



^/l _ e -2/f(t-ii) 



+ R N+2 (t,y). 



By (3.31), we deduce that 
^lC TJ B|/C(t,xo,y) 

TV 



Y^{yKAC T B^ 2 Vi) )AB 3 T /2 Vfrre KT )C(t,xo,y) jf _ === dU 



+ R N+2 (t,y). 



Following (3.33), we have: 

(is 



vT 



-2K(t-s) 



V2KT 



1 - exp(-2KT) 7 v 72 ^^ - s 



T 



'1 - exp(-2KT) Jo y/t-s 



'1 - exp(-2ATT) 



Since CV = v 27r exp(f^T)CV, we deduce that (3.35) holds at rank N + 1. 
Third Step. Suppose that VI < [v^ACr^ 2 ]' 1 / 2 - Then > the 

series in (|3.35l) is 

convergent (and bounded above by 2). The point is thus to prove that -Rtv(£, y) — > 
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as iV — > 00 for all y ^ 1. By Lemma 3.6, we deduce that, for \/t ^ L4ki(T)] 1 



RN(t,y) 



N-l 



o < t N «: «: ... ^ t v^v v 7 ! - e 2E "(^ 



zfitx ... (it 



A' • 



Using (3.26), we obtain 



RN(t,y) ^ y2ixAB, 



,3/2 



vT 



Ak\{T) 



N-l 



x / — — JJ 1 dti 



AT- 



(3.36) 



To compute the integral in the right-hand side of the above, define 



J n (t) = I u {n - 1)/2 {t - u)- 1/2 du, £ >0, neN. 
Jo 



By the change of variable u = ts, we have J n (t) = t n / 2 J n {l). Note that J$(l) 



7r ^ 4. Moreover, for n ^ 1, 



1-1/n 



n + 1 



(3.37) 



With this notation 



AT-1 



r= TI /f .\. d *i • • • 



< tjv < tjv-i < . < to V ^JV i=Q ■ H+l 

N-2 



< tjv-i < ... < to 



Jo(tN-i) TT ,. , =dti . . . dt N _i, 



i=0 



where the second equality holds for N ^ 2. More generally, setting for any n ^ N—2 



lN,n(t) :- 



< tjV-(n + l) < ■■• < to 



N-(n+2) 

Jn(tN-{n+l)) I 1 , =^1 • • • dtpf— 

V ^ — 



+1)' 



i=0 
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we have, for any ^ n ^ iV — 3, 



I N ,n{t) = J«(l) 



t 



< i]v_( n+ i) < ... < to 



n/2 

iV-(n+l) 



jV-(ra+2) 

n 



i=0 



i+l 



-.dt\ . . . dt^_^ n+ i) 



'0 < ijv_( n+ 2) ^ ••• ^ *0 

7V-(n+3) 



JV-(n + 2) 



n/2 

"iV-(n+l) 



y/tN-{n+2) — £/V-(ra+l) 



AT-(n+l) 



n 



i 



i=0 

<^n(l)-^A r ,n+l(^), 



i+l 



. . . dtN-(n+2) 



so that, by (3.37), 



riV-3 



n 



n=0 
■JV-1 



{ N,N-2 



(t) 



•N-2 



n 



n=0 



JiV-l(t) 



n 



L n=0 



using the fact that 7^,^-2 (*) = «/jv— 2 (1) «7?sr— i (^) • Going back to (3.36), we deduce 
that R N (t,y) — > as N — > oo. Hence taking the limit as iV — >■ oo in (3.35) yields 
(with y/t < [v / i?lC i TJ By /2 ]- 1 /2) 

|i^(t,y)| ^2lC r flfX(f,x ,!/), 
from which the result follows. □ 

Complementary to the above two results, we have the following two, which esti- 
mate the difference between the densities of two processes driven by two different 
drifts. 

Lemma 3.8. Let T > 0. Suppose Xq — xq < 1, and that ei,e2 G £(T,A). Then 
there exists a constant Ks(T) (independent of A and Xq) which increases with T such 
that 

sup [yG\\d y pl°(s,.) -d yP x e °(s, 0||J ^^c 3 (T)||ei-e , 2 || 00>t , 

< s < * 

for allt^ min{[2K 3 (T)]- 2 ,r} ; where A := max(aA + \b(l)\,l). 



Proof. Let y ^ 1. Again by (3.10), we see that 



<(aA+\b(l)\) 
+ «l|ei -e'alloct 



./-oo 



|&p*°(s,;z) -9*p^(s,s)| |0 B g(t-s,s,y)|tfeds 



d z pl°(s, z) - d z q(s,x ,z)\ \d y q(t - s, z,y)\dzds 



J-oc 

t r ± 



II 1 ziioo.i 



|<9 2 g(s, x , z) | |<9 y g(i - s, z, y) \ dzds. 



(3.38) 



J-oc 
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Using Proposition 3.5 once again and the notation (3.30) (see also (3.28) and (3.29)) 
together with (3.26), we see that 

lW.»)-W,y)| 

*C Ak 3 (T) sup [v^ \\d z p*°(s, .) - d zP :°(s, .)||J 

\d z Peo( s , z ) ~ d z q(s,x ,z)\ 



sg s sg t 



t rl 



+ k 3 (T) \\e[ - 411^ 
+ « 3 (T) llei-eall^^ y 



J-oo 



/C(— (i — s), 2, y)dzds 



1 1 



=/C(s, xq, z)K{— (t — s),z, y)dzds, 



fs y/t — s 

for some constant k 3 (T) > 0, where A = max(aA + |6(1)|,1). Using Proposition 



3.7 and (3.31), and allowing the constant K 3 (T) to increase as necessary from line 



to line below, it follows that 
\d y p%%y) - d y pH{t,y)\ 

^Ak 3 (T) snp [^||«( S ,.)-«(s,-)||J 



< s < i 



+ K3 (T)A\\e' 1 -e' 2 



2 lloo,f 



J-oo 



/C(s, x , 2)/C(— (£ — s), z, y)dzds 



+ k 3 (T) \\e[ - e'all^ t /C(t, x , J/) 



for all t ^ min-fLAfi^^)] 2 ,T}. Using (3.31) again, together with the bound A ^ 1, 
we deduce: 

\d y p^(t,y)-d yP T 2 (t,y)\^AK 3 (T) sup [fi \\d y p x e °(s, .) - d y p%(s, .)|| J 

< s < t 

+ k 3 (T)2 He; - 411^ £(t, x , j/). (3.39) 

Assuming that K 3 (T) ^ K^CO and multiplying the above inequality by y/i, we 
deduce that, for t ^ min{[2AK3(T)]~ 2 , T} we have 

sup [^\\d y p:°(s,.)-d y p*°(s,.)\\j 



^ s ^ t 



^ 2k 3 (T) sup 

0<t ^ T 



B T (exp(2Kt) 



AW' I 

A \\e 1 — e 2 



2 lloo.t ' 



□ 



Proposition 3.9. Let T > 0. Suppose Xq = xq < 1, and that ei,e2 € C(T,A). 
Then there exists a constant K^iT) (independent of A and xq) which increases with 
T such that 

|2 S 



1 



\d y p%%y) - d y p2(t,y)\ < Ak a (T) — exp 



lff°-y| 

K A {T)t 



'1 e 2lloo,t) 



/or all t ^ minjfAK^T)] 2 ,T} and y ^ 1, where (£,t°)t>o ^ s as ^ n Proposition 
and A := max(aA + |6(1)|, 1). 



3.5 
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Proof. We follow the strategy of the proof of Proposition 3.7| We thus define 

GT 1 , e2 (t,y) = d yP T 1 (t,y)-d y pT 2 (t,y). 



Going back to (3.38) we can proceed in the same way as in Lemma 3.8 in order to 
bound the second two terms of this expression. We then get that, for some constant 
ka(T) > 0, 



fe 2 (ti,yi)| 



+ Ak a {T) 



y/\ - e -2A-(t-ti) 



K.(-(t - t 1 ),y 1 ,y)dy 1 dt 1 . 



We can iterate this inequality in exactly the same way as in Proposition 3.7 (pre- 
cisely, we can divide both sides by ||ei _ e^lloot * n OT< ^ eT t° recover (3.34) with 
\F*°(-,y)\ replaced by \G^° e2 (-, y)\/\\ei — e' 2 \\oo,t therein and then follow (3.35)), and 
the convergence follows from Lemma 3.8 This yields the result. □ 



3.3. Conclusion. 



Proposition 3.10. Let T > and X be as in the statement of Theorem 3.1 with 
(w.l.o.g.) (3^1. Suppose e^e-i 6 T-L(T,A,Xq). Then there exists a constant 
k(T), independent of A, (5 and e, and increasing in T, and a constant K,(T,/3,e), 
independent of A and increasing in T , such that 



sup 

< s < t 



< AK(T,(3,e)Vi(l +E\X \ 



fort < mm{[AK(T)}- 2 ,T}, where A := max(aA + |6(1)|,1) 



Proof. We have by fl2.4|) and (|3.6|) 

d 



dt 



E(M t ei ; 



d 
dt 



E(Mf 



1 



\d y pl(t,l)-d yP : 2 (t,l)\F(X edx) 
\ E f ^vP^it " M) - dyP^it ~ s > 1)1 P K ei G da) 

(Li + L 2 + L 3 ). 



(3.40) 



3.9 



and 



Suppose t ^ T and y/l ^ [A/? 4 (T)] 1 , where K&(T) is as in Proposition 
A = max(aA + |6(1)|, 1). Considering the first term only, we can use this result to 
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see that 

L x ^ Af3K A (T) ( 

+ Ak 4 (T) 



\Ji-e y/t 



exp 



(1-x) 2 

K A (T)t 



(1 - x)dx J || e\ - e' 2 \\oo,t 



l-e 



Vi eW { K,{T)t) 



F(X e dx) )\\e' l -e 2 



I oo,i j 



where we have used the bi-Lipschitz property (|3.19|) in the statement of Proposi- 
tion 



3.9 (and thus modified the value of n<i(T) accordingly), together with the fact 
that £l = 1 for all t ^ 0. Indeed, by definition (see the statement of Proposition 



3.5), we have 



KQds, 6(1) = 0. 



We deduce that there exists a constant /ts(T, 0, e) > 0, which is independent of A 
and which is allowed to increase as necessary from line to line below, such that (the 
value of k^T) being also allowed to increase from line to line) 



Li ^ ApK 4 (T)Vi ( / zexp 



dz ) || e[ — e' 2 



2 \\oo,t 



+ Ak a (T) — exp 

y/t V K 4(J )t 

^AK S (T,0,e)Vt\\ 



'1 e 2ll°o,t 



(3.41) 



We can then use Proposition |3.9| again to see that 

1 



L 2 ^ Ak 4 (T) sup 

0<s ^ t 



s ^ 2 exp 



K 4 (T)s 



E(M?)\\e[-e', 



2\\oo,t- 



By Proposition |2.2[ we deduce that 

L 2 ^ A.K A (T)y/i(y/i + E\X \) \\e[ - e' 2 



2l|oo,t) 



(3.42) 



where we have used the elementary inequality exp(— l/v) ^ v for all v ^ 0. We 
finally turn to L3 in (3.40). By Proposition 3.7, Proposition 3.5 and (3.19), we also 
have that 



1 



\d v p%'\t-s,l)\ < K 4 (T)A-==exp 
^ Ki(T) exp 



1 



K 4 (T)(t-s) 
A 



n 4 (T)(t-s) 



\d y q(t-s,0,l)\ 
1 



y/t - S t — S 



< k 4 (T)A. 



(3.43) 
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Thus, from fl3~40j ), fl3~4lj ), fl3~lL2| and fl3~43] ), we deduce 

^ lfi: 5 (T,/3,e)v / t(l + E|Xo|)||e / 1 -e'^l^ 
d 



e 2 ■> 



ds. 



By taking the supremum over all s ^ t in the above, we have, for t ^ (2k 4 (T)/1) 1 
(which actually follows from the aforementioned condition t ^ (k 4 (T)v4)~ 2 by as- 
suming w.l.o.g. /t 4 (T) ^ 2), 



sup 

o c s < t 



^ 2Ak 5 (T, 0, e)Vt (1 + E|X |) He; - e' 2 \\^ t . 



□ 



Proof. [Proof of Theorem 3.1 Choose 

Ai = 2 sup 

o ^ t ^ i 



E(M° 



and set 

A 1 = max(aA 1 + |6(1)|, 1). 
Note that Ai and Ai depend on j3. Then choose T x ^ min {[^(l)]" 2 , 1} such that 

T\^(l,/3,e)l 1 (l + E|X |) ^ ], (3.44) 



where k(1) and ft(l,/3,e) are as in the Proposition 3.10 By that result, if e G 
ft(A,7i,X ) then 



E(M t e 



*C >/tK(Ti, /3, e)A(l + EIXqDA + 



E(M t ° 



for all t ^ min{[74i/t(Ti)] 2 , 1}. By definition, we have T\ ^ 1 so that k(Ti) ^ k(1) 
and k(Ti, (3, e) ^ «(1, /3, e). Therefore 



E(M t e 



^ VtK(l,P,e)A 1 (l + E\X \)A 1 



E(M t °; 



for all t ^ min{[Ai/t(l)]- 2 , 1}. Hence for all t G [0,Ti 
d 



E(M t e 



^ >/t«(l, /3, e)A x (l + E\X \)A 1 + 
d 



dt 



E(M t °) 



< — + sup 

* < i < 1 



df 



e(m°; 



by ggp , so that T(e) G ^(A, Ti, X ). 

To prove that T is a contraction on TL(Ax, T 1? X ), first note that for e G "H(^4i, Ti, Xq) 

||e'||oo,Ti ^ ||e||'H(A,ri,Xo) ^ 2||e'|| 00i r 1 
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by the mean value theorem, since e(0) = and T\ ^ 1. Thus for any ei,e2 G 
H(Ai,T x ,X ) 

\\Y{e x ) - T(e 2 )\\ n{AuTuXo) <: 2 ||r( ei )' - IWU^ 



^ 2 sup 



If (MP) - ±MM, 



^ 2^/T\K{T 1 ,f3,e){l+E\X \)A 1 \\e , 1 -e , 2 \\ OO;Tl 
1 

< d| e i - e2\\H(Ai,Ti,x ), 



by our choice of T\ and using Proposition 3.10 once more. Since T-L(Ai, T\, Xq) is a 



closed subspace of C 1 [0,T] (a complete metric space), the existence of a fixed point 
for T follows from the Banach Fixed Point Theorem. □ 



4. Long-time estimates 

In order to extend the existence and uniqueness from small time to any arbitrarily 
prescribed interval, we need an a priori bound for the Lipschitz constant of t i— >■ 
E(M t ) on any finite interval [0, T\. In the whole section, for a given initial condition 
Xq = xq G (— oo, 1), we thus assume that there exists a solution to (2.1) according 
to Definition 2.3 i.e. such that e : [0, T] 3 1 1— >■ E(M t ) is continuously differentiable. 
The aim is to then prove a universal bound for the Lipschitz constant of [0,T] 3 
t h-> E(M t ): 

Theorem 4.1. For a given e G (0,1), there exists a positive constant a® G (0,1], 
only depending upon e, the Lipschitz constant K and the linear growth constant A 
of b, such that, for any a G (0, cto) and any positive time T > 0, there exist a 
constant M.t, only depending on T , e, K and A, such that, for any initial condition 
Xq = Xq ^ 1 — e, any solution to (2.1) according to Definition 2.3 satisfies 

WG[0,T], e'(t) = jE(M t )^M T . 

Remark 4.2. In the proof below, we provide an expression for oq (or put it differ- 
ently a lower bound for the best possible choice of ao): we refer to the statement of 
Proposition ^. 5 , where the constant cxq therein is the same as the one in the statement 
4-l\ Without any further assumptions on the coefficients, the expression 



of Theorem 
in Proposition 

stant. However, it becomes tractable when b(x) = —Ax, which is the common case 



4-5 is not completely explicit as it depends upon a non-explicit con- 



treated in the neuroscience literature. We thus refer the reader to Subsection 6.1 for 
a detailed discussion about the critical threshold ao in this framework. 

4.1. Reformulation of the equation and a priori bounds for the solution. 

In the whole proof, we shall use a reformulated version of the nonlinear equation 



(2.1). Indeed, given a solution (X t , M t )o <g t ^ t to (2.1) on some interval [0, T] ac- 



cording to Definition 2.3, we set 

Z t = X t + M t , te[o,T). 



(4.1) 
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Then, the dynamics of (Z t )o ^t^T have the form: 



b(X s )ds + aE(M t ) + W t , te [0,7], 



(4.2) 



so that, by continuity of the mapping [0,7] 3 t i— >■ e(t) = E(M 4 ), the process 
{Z t )o^t^T has continuous paths. Moreover, as noted above in Proposition 2.2 Mt 
is completely defined in terms of Z t , since it holds that 



(4.3) 



Lemma 4.3. There exists a constant B(T,a,b), only depending upon 7, a, b and 
non- decreasing in a, such that 



M t =[( sup Z s )\ = sup [{Z s ) + \. 
Following the proof of Proposition |2.2[ we claim: 



sup e(t) = e(T) ^ E 

sS t ^ T 



sup (Z t ) + 

t sg T 



^B(7,a,b). 



(4.4) 



A possible choice for B is 



r>(rp t\ E[(X ) + ] + 4T 1 / 2 + AT 
5(T, a, b) = exp 



a: 



2AT 



a: 



Proof. The proof is a straightforward consequence of (2.12) in the proof of Proposi- 
tion [221 □ 

4.2. Local Holder bound of the solution. We now turn to the critical point of 
the proof. Indeed, in the next subsection, we shall prove that, for a small enough, 
the function t i— > e(t) = E(M t ) generated by some solution to (2.1) according to 
Definition 2.3 (so that e is continuously differentiable) satisfies an a priori 1/2- 
Holder bound, with an explicit Holder constant. This acts as the keystone of the 
argument to extend the local existence and uniqueness result into a global one. As 
a first step, the proof consists in establishing a local Holder bound for e in the case 
when the probability that the process X lies in the neighbourhood of 1 is not too 
large. 

Lemma 4.4. Keep the same notation as in the previous subsection and consider in 
particular a solution (X t )o^t^T to (2.1) on some interval [0,7], with 7 > and 
initial condition Xq = xq. Assume in addition that there exist some time to G [0,7] 
and two constants e G (0, 1) and c G (0, 1/a) such that Xq ^ 1 — e, and, for any 
Borel subset Ac [1 — e, 1], 

F{X t0 G A) ^c\A\, (4.5) 
where \A\ stands for the Lebesgue measure of A. Then, with 

exp(2A) [(8 + 5c + 8e _1 )A + 4(2 + c + e" 1 )] 



1 — ca 



it holds that, for any h G (0, 1), 
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Proof. Without any loss of generality, we can assume to = 0, with T being under- 
stood as T — to. Indeed, setting 

Xl t0 := X t0+t , te[0,T-t o ], (4.6) 
we observe that, for t G [0, T — £ ]j 

xl t0 =X t0 + b{X$ t0 )dr + aE{M t+t0 -M tQ )+W t+ t -W to -(M t+t0 -M t0 ). (4.7) 
Jo 

Here M t+to — M to represents the number of times the process X reaches 1 within the 
interval (t ,t+to}. Therefore, this also matches the number of times the process X"'o 
hits 1 within the interval (0, t], so that X^o indeed satisfies the nonlinear equation 

(2.1) on [0, T — to], with X '° = X to as initial condition and with respect to the 

shifted Brownian motion (W t tQ := W to +t — Wto)o<t <T-t - m wna t follows, t is 
thus assumed to be zero, the new T standing for the previous T — to and the new 
Xo matching the previous X to and thus satisfying (4.5). 

For a given h G (0, 1), such that h ^ T, and a given 23o > (the value of which 
will be fixed later), we then define the deterministic hitting time: 

R = 'm£{t G [0, h] : E(M t ) = e(t) ^ B h 1/2 }. 



Following the proof of (2.11) (see more specifically (2.9)), we have, for any t G 
[0,hAR], 

M t < sup (Z s ) + (X ) + + A / (1 + (Z s ) + + M s )ds + ae(t) + 2 sup 

Osgssgt Jo sg s sg t 

^(X )+ + 2a/ (l + M s )ds + aS /i 1/2 + 2 sup |W S | 

JO < s < i 

^ (X )+ + 2A/i + 2A / M s ds + aB h 1/2 + 2 sup |W,|, 

JO ^ s < t 

where we have used (4.3) to pass from the first to the second line. By Gronwall's 
Lemma, we obtain 



M t ^ exp(2A/i) 

< (X )+ + exp(2A/i) 



(X )+ + 2Ah + aB h 1/2 + 2 sup |W, 

^ s ^ h 



Akh + aBoh 112 + 2 sup |W S 

^ s^h 



(4.8) 



as exp(2A/i) < 1 + 2A/iexp(2Ah) and (X )+ < 1. 

Assume that Bq exp(2A/i)/i 1//2 ^ e/2 ^ 1/2. Then, by Doob's L 2 inequality for 
martingales, 



P(M t ^ fc) < P ( exp(2A/i) 



k > 2 



fc > 2 



4A/i + 2 sup |W S 



^k-3/2 



(4.9) 



< 2exp(2A/i)E[4A/i + 4|W4|] < exp{2Ah) [8Ah + 8/i 1/2 ] . 
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Moreover, 
P(M t ^ 1) 



^ P I (AT )+ + exp(2A/i) 
^ P ( X e [1 - e, 1],X Q + exp(2A/i) 



AAh + aB h 1/2 + 2 sup \W S 

^ s sg /i 



^ 1 



4A/i + ai3 /i 1/2 + 2 sup |VT S 



> 1 



AAh + 2 sup |W S 

s£ s s£ fi 



^ e/2 



+ P ^exp(2A/i) 
:=/i + / 2 , 

where we have used i3 exp(2A/?,)/i 1 / 2 ^ e/2 in the third line. 
By Doob's maximal inequality, we deduce that 

J 2 2exp(2A/i)e~ 1 E[4A/i + 2|W ?t |] exp(2A/i)e~ 1 [8A/i + Ah l/2 ] . (4.10) 

We now switch to Jj. By independence of X and (Ws)o <: s < r and by (4.5), 



Ji < c / P exp(2A/i) 



4A/i + a£ /i 1/2 + 2 sup |W, 





^ C 



P ( exp(2A/i) 



4A/i + a£ /i 1/2 + 2 sup |W S 



^ xj <ix 
^ x ) dx 



AAh + aB h 1/2 + 2 sup |W S 

o s; s s£ h 



= cexp(2A/i)E 

By Doob's L 2 inequality, 

h <: cexp(2A/i) [4A/i + aB h 1/2 + Ah 1/2 ] . 
Together with ( |4.10[ ), we deduce that 

P(M t ^ 1) < exp(2A/i) [4(c + 2e~ 1 )A/i + 4(c + e" 1 )^ 1 / 2 + caB h 1/2 ] . 
From ( 4.9[ ), we finally obtain, for t ^ R A h, 

E(M t ) = W > k ) 

^ exp(2A/i) [4(2 + c + 2e~ 1 )A/i + 4(2 + c + e^ 1 )/^ 2 + caB h 1/2 ] 

<: exp(2A/i) [(8 + 5c + Se'^Ah + 4(2 + c + e^ 1 )/^ 2 ] + caB Q h 1 ' 2 , 

provided B$ exp(2A/i)/i 1 / 2 ^ e/2 ^ 1/2, which implies 

ca£ exp(2A/i)/i 1/2 ^ caB h 1/2 + cA/i, 

using the fact that exp(2A/i) ^ 1 + 2Ah exp(2A/i). Therefore, if R ^ h, then we can 
choose t = .R in the left-hand side above. By continuity of e on [0, T], it then holds 
e(R) = Boh 1 / 2 , so that 

(1 - ccO£ /i 1/2 < exp(2A/i) [(8 + 5c + Se'^Ah + 4(2 + c + e^h 1 ' 2 ] 

< exp(2A) [(8 + 5c + 8e _1 )A + 4(2 + c + e" 1 )] h 1/2 , 
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which is not possible when 

exp(2A) [(8 + 5c + 8e _1 )A + 4(2 + c + e" 1 )] 



ca 



Precisely, with Bq as above and Bq exp(2Ah)h 1 ^ 2 ^ e/2 it cannot hold R ^ h. 



4.3. Global Holder bound. In this subsection, we shall prove: 



□ 



Proposition 4.5. For a given e G (0, 1), there exists a positive constant Oq G (0, 1], 
only depending upon e, K and A, such that, for any a < «o, there exists a constant 
B, only depending on a, e, K and A ; such that, for any positive time T > and any 
initial condition Xq = xq ^ 1 — t, any solution to (2.1) according to Definition 2.3, 
so that the mapping e : [0, T] 3 t i— >■ E(M t ) is continuously differentiable, satisfies, 
for any h G (0, 1) and to G [0, T], 

Aote £/iat i3 above may differ from Bo in the statement of Lemma 4-4 The constant 
ciq can 6e described as follows. Defining T as the largest time less than 1 such that 

(l - e) exp(ATo) ^ 1 - 7e/8, AT exp(AT ) < e/8, 

then ao can be chosen as the largest (positive) real satisfying (with B(To,ao,b) as in 
Lemma 4-3) 



a B(T ,a ,b) ^ e/4, 
3/2 — (~) 



ao2 3/2 (c'j"'-exp 



B{T ,a ,b)} < 1, 



C i 



2 3 / 2 ( c '; 



3/2 



exp 



B(T ,a ,b) 



^ 1, 



where the constant d is defined by the following property: d > depending on K 
only is such that for any diffusion process {U t )o ^t^i satisfying 

dU t = F(t,Ut)dt + dW t , £G[0, 1], 

where Uq = and F : [0, T] x R — y M is K -Lipschitz in x and satisfies F(t, 0) = 
for any t G [0, 1], it holds that 



d d 
—nU t G dx) ^ —= exp 

dx y/t 



X 

di 



xeR, te (0,1]. 



Xq ^ 1 — e, with e G (0, 1), and a 



The proof relies on the following 

Lemma 4.6. Given an initial condition Xq 
solution {X t )o^t^T to (2.1) on some interval [0, T], the random variable X t has a 
density on (— oo, 1], for any t G (0,T]. Moreover, defining T as in the statement of 
Proposition 4-5\ and choosing a ^ a\ satisfying 

ai5(T ,ai,6) ^ e/4, 
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then, for x£ [1 — e/4, 1), 



d 

dx 
d 

dx 



P(X t GdxK2 3 /V) 3/2 exp(-0 [e- 1 + B{T , a, &)] ift^T 
F{X t G dx) <: cX 1/2 + 2 3/2 (c) 3/2 exp (-^j B(T ,a,b) if t > T 



where the constant d is as in the statement of Proposition 4-5 



Before we prove Lemma |4.6[ we introduce some materials. As usual, we set 
e(t) = E(Mf), for t G [0, T], the mapping e being assumed to be continuously 
different iable on [0,T]. Moreover, with pQ)o<*<T) we associate the sequence of 
hitting times (rjj,)^ o given by (2.3 ) . We then investigate the marginal distributions 
of (Xt)o <t <r- Given a Borel subset A C (— oo, 1], we write in the same way as in 



the proof of (3.1 ) 



F(X t eA) = F(X t eA,n> t) 



k >i J ° 



(4.11) 



where the notation P(-|rfc = s) stands for the conditional law given = s. Following 
( 4.6[ ) and (4.7), we can shift the system by length s G [0,T]. Precisely, we know 
that (X| s := X s+r ) ^ r ^ t- s satisfies 

pr 

XI s =X S + b(Xt)du + ae is (r) + W s+r - W s - M» s , (4.12) 
Jo 



with 



e «-(r):=e(s + r)-e(s) ; M}' := ^ l (0 , r] (t£) ; 

k>i (4.13) 

rt :=M{u>Tt 1 :X a+Ur _>l}, k^l ; r«* := 0. 

Conditionally on = s, the law of (X| s ) o ^ r < T-s until rf s coincides with the law 
of (Z^ s,0 ) ^ r ^ y-s until the first time it reaches 1, where, for a given J^-measurable 
initial condition £ with values in (— oo, 1), (Z^. s ^) ^ r ^ y_ s stands for the solution of 
the SDE: 



Z^^ = ( + b{Zl s ^)du + ae h (r) + W r , re[0,T-s}. 



(4.14) 
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Below, we will write Z£ for Z^ ^. In the end, we thus have (by (4.11)): 
F(X t eA)=F(X t eA,r l >t) 

[ HXts e A rt >t-s\r k = s)F(r k G ds) 



k > 1 



X 



■4) + E/ 



P(4-° G A)P(r fe G ds) 



(4.15) 



fcjs 1' 
t 



= F(Zf° G A) + / P(Zh u € A)e'( S )d S , 
Jo 

for any Borel set A C (— oo, 1], the passage from the third to the fourth line following 
from (3.9). 

We can now turn to: 

Proof of Lemma \4.6] Given an initial condition Xo G (— oo, 1 — e] for e G (0, 1), we 
know from Delarue and Menozzi [6j that Z^° has a density for any t G (0, T] (and 
thus Zfl!g as well for ^ s < t). From (4.15), we deduce that the law of X t has 
a density on (— oo, 1] since F(X t G A) = when \A\ = 0, where \A\ stands for the 
Lebesgue measure of A. Moreover, there exists a constant d ^ 1, depending on .fT 
only, such that, for any t G [0, T A 1]: 

d - c' 
-F(Zf° G dx) ^ — exp 



dx 



where is the solution of the ODE: 



\x 



dt 



d 
dt 



& t = b(0 t )+ae?(t), t G [0,T], 



(4.16) 



(4.17) 



with = x . Above, the function [0, T] 9 t H> e(t) represents [0,T] 9 t (-)> E(M t ) 
given Xo = Xo, which means that the initial condition Xo of Xo upon which e depends 
is fixed once and for all, independently of the initial condition of i?. In particular, 
as the initial condition of $ varies, the function e does not. We emphasize that d 
is independent of e and can be taken to be that defined in Proposition 4.5 Indeed, 
we can write P(£f° G dx) as P(Z£° - 0*° G d(x 



0f )), with 



d(Z?-0?>)=F(t,Z : 



■d x t °)dt + dW u te[0,T\ 



ryXQ 



a 



F(t,x) = b(x + ti%°) tG[0,T], xG 



We then notice that F(t, •) is X-Lipschitz continuous (since b is) and satisfies 
F(t, 0) = 0, so that, referring to |6], all the parameters involved in the definition of 
the constant d are independent of e. The fact that d is independent of e is crucial. 
As a consequence, we can bound (d/dx)F(Zl s 2s G dx) in a similar way, that is, with 
the same constant d as in (4.16): for any ^ s < t ^ T, with £ — s ^ 1, 

d 



dx 



P(i{i{ G dx) < 



tU,0l2 



exp 



c'(t 



(4.18) 
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where & s, ° is the solution of the ODE: 



d f=b(4°)+a^(t), te[0,T-s) 



with 



dt 

as initial condition. 



Bound of the density in small time. Keep in mind that X = x ^ 1 — e. Therefore, 
by the comparison principle for ODEs, ^ for any t G [0,T], so that by 
Gronwall's Lemma 

d x t ° ^ d]- e < (l - e + AT + ae(T)) exp(AT). 



By Lemma 4.3, we know that e(T) ^ B(T,a,b), so that 

tff° ^ (1 - e + AT + aB(T, a, b)) exp(AT). 



(4.19) 



Now choose T as in Proposition 4.5, i.e. T ^ 1 such that 

(1 - e) exp(ATo) ^ 1 - 7e/8, AT exp(AT ) ^ e/8, 
and then take ot\ G (0, 1) such that 

ai B{T Q , a u b) exp(ATo) < e/4. 
Then, whenever a. ^ oti, it holds that 

t G [0, T A T]. 

Therefore, for a; ^ 1 — e/4, 



cxp 



^012 



exp 



16c't 



t G [0, T AT]. 



(4.20) 



Similarly, 



i?K ^ 3e/8 ^ 3/8, 0O^^ AT. 



Indeed, e^(T-s) ^ e(T) for s G [0, T], so that Q4.19| ) applies to t?{fj° with 1-e therein 

tt 

being replaced by 0. Therefore, for x ^ 1 — e/4, it holds that x — w^lg ^ 3/4 — 3/8 = 
3/8 ^ 1/4, so that 



tt s ,o l2 



cxp 



c'(£ - s) 



exp 



16c'(t - s) 



^ s < t C T n A T. 



(4.21) 



In the e nd, fo r x G (1 - e/4, 1) and t < T AT, we deduce from ( |4.15[ ), ( |4.16D , ( |4.18[ ), 
(4.20), (4.21) and Lemma 4.3 again, that 



d 

dx 



¥{X t G dx) ^ c'ujq [e- 1 + e(T A T )] ^ c'w [e' 1 + B(T Q , a, b)] , (4.22) 



where 



wq = sup 

i>0 



t~ 1/2 exp 



Wc't 



4v / c 7 sup[wexp(-w 2 )] = 2 3/2 ^c'exp ( -\ ) . 
«>o \ 2/ 
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Bound of the density in long time. We now discuss what happens for T > Tq and 
te[T ,T}. Then, 

£ppr, 6 «fa) < 6 rf'- < r„) + A P ( Xl 6 ^ Tl »- > To ) (4 23) 

= 7Ti + 7T 2 , 

with rJ*- T ° = inf{w > : X t _ To+u _ ^ 1} = inf{w > : X^Z T ° > 1}. The above 
expression says that we split the event (X t is in the neighbourhood of x) into two 
disjoint parts according to the fact that X reaches the threshold or not within the 
time window [t — Tq, t]. We have chosen this interval to be of length To in order to 
apply the results in small time. 

We first investigate 7r 2 . The point is that, on the event that rf t_T ° > T and within 
the time window [t — T ,t], X behaves as a standard diffusion process without any 



jumps, namely as a process with the same dynamics as Z^- T o' Xt - T o , Following (4.16) 
we then have 



< ip (z^ o - To ' Xt - To G dx) < sup ip (4r T °' 2 e dx) < c'T 1/2 . (4.24) 



We now turn to 7T\. Here we write that 



7T! = £P (X t G dx, rf^ < T ) = £ J-P G ^ r »-o ^ Tq < r «^,) 

fe > 1 

= E 1 P ( X * e ^~ T ° < T ° < \-t T ° = s ) nrt To e ds) 

= E f 1 F (^-; T °'° e dx,r < ft?) Hrt T ° G ds), 
fc > 1 "' 



since on the event that r fc f T ° ^ T < t^,^° and given that the fc-th (and last) jump 
of X in the interval [t — T ,t] occurs at time t — T + s with s G [0,T ], we have 

that the process X r for r G [t — T + s,t] coincides with the process zt +t ~ T °'° for 
w G [0, T — s]. Thus 



*i< E 

fc> i 



t ° d 

dx 



f ( G dx) P(^^ G d.s) 



(4.25) 



o 



— P [Z r 

dx V r -s 



tts+t-T .0 



G dx ) e'(s + £ — T )ds. 
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By ( |4.18[ ), we have 
rTo d 



d /(z- + _!; T °'"edx)e'(s + t 



To 



exp 




e'(s + 1 — T )ds. 



Following the shift in (4.12), it is well seen that the mapping [0,T ] 9s4 e^~ T o(s) 
satisfies Lemma [473 , that is 



sup e* 

o ^ s ^ T 



r o( s )= sup [e(s + t-T )-e(t-T )]=e(t)-e(t-T )^B(T ,a,b). 

s: s sg T 



Therefore, we can follow the same strategy as in short time, see (4.21) and (4.22). 
Indeed, for a ^ a%, by the choice of To as before, it holds that 

7Ti ^ c'w B(T ,a,b), 



for x G [1— e/4, 1). Using (4.24) and the above bound, we deduce that, for t G [To, T], 



d 

dx 



P (X t G dx) < c' [T " 1/2 + ^o5(T , a, 6)] • 



□ 



Now, we can complete the proof of Proposition 4.5 



Proof of Proposition J^5_. Proposition 45 follows from the combination of Lemmas 
and 46 Indeed, given To and «o as defined in Proposition 4.5[ then by Lemma 
it follows that ¥(X t G A) < (l/a)\A\ for any Borel subset A C [1 - e/4, 1], any 



4.4 



4.6 



a < «o and any t G [0, T]. The result follows by Lemma 4.4, with B being given by 
i3oexp(2A) with e in Bq replaced by e/4. □ 

4.4. Estimate of the density of the killed process. In light of the previous 
subsection, for a solution (X t )o to (2.1) such that the mapping [0,T] 3 t i-> 

e(t) = E(Mj) is continuously differentiable, we here investigate 

— P (X t G dx, t < n) , t€[0,Tl, x^l, 

where 7"i = inf{t > : X t _ ^ 1} as usual. This is the density of the killed process 
(X tAT1 ) <^ t ^T, which makes sense because of Lemmas 3.3 and 4.6 
Here is the main result of this subsection: 

Lemma 4.7. For any e G (0, 1), any T > and any B > 0, there exist two constants 
Ht andrjT, only depending upon T , B, e, K and A, such that, for any initial condition 
Xq ^ 1 — e and any continuously differentiable non-decreasing deterministic mapping 
[0, T] 3 t i-> e(t) satisfying 

e(0) = ; e(t) - e(s) < B(t-s) 1/2 , O^s^t^T, 

if (xt)o < t < t denotes the solution of the SDE 

dXt = b{xt)dt + ae'(t)dt + dW t , t G [0, T] ; xo = £o, 
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then, with p(t,y) denoting the density of Xt killed at 1 as in (3.3), we have 

p(t,y)^fi T (l~yr, t G [0,T], ye [l-e/4,1]. (4.26) 

Proof. First Step. The first step is to provide a probabilistic representation for p. 
For a given (T,x) G (0, +oo) x (—00, 1), we consider the solution to the SDE: 



dY t = - [b{Y t ) + ae'(T - t)] dt + dW u t G [0, T] ; Y = y, 



(4.27) 



together with some stopping time p ^ p A T, where p = inf{t G [0, T] : Y t ^ 1} 
(with inf = +00). Then, by Lemma 3.3 and the Ito-Krylov formula (see [T2l 
Chapter II, Section 10]), 

d(p(T-t,Y t )) 

= -d tP (T - t, Y t )dt - [b(Y t ) + ae'(T - t)] d y p(T - t, Y t )dt + \d 2 yy p{T - t, Y t )dt 

+ d y p(T-t,Y t )dW t 
= b'(Y t )p(T - t, Y t )dt + d y p{T - t, Y t )dW t , 
for ^ t ^ p. Therefore, the Feynman-Kac formula yields 

(4.28) 



p(T,y)=E p(T - p,Y p )l {Ypm exp(- J\'(Y s )ds\ , 
the indicator function following from the Dirichlet boundary condition satisfied by 

Second Step. We now specify the choice of p. Given some free parameters L ^ 1 
and 5 G (0, e/4) such that L5 ^ e/4, we assume that the initial condition y in (4.27) 
is in (1 — 5, 1) and then consider the stopping time 

p = M{t G [0, T] : Y t <£ (1 - L5, 1)} A 5 2 . 



Assume that 5 2 ^ T. By ( |4.28[ ), we deduce that 



p(T,y) ^ exp(/i-5 2 )(l-P(y p = l)) sup p(t,z), 

(t,z)eQ{S,L) 



(4.29) 
(4.30) 



with 



Q(5, L) = { (t, z) G [T - 5 2 , T] x [1 - L5, 1] } . 

The point is then to give a lower bound for F(Y p = 1). By assumption, we know 
that e is i3-H6lder continuous on [0, T]. Therefore, since Y — y G (1 — 5,1), we 
have, for any t G [0,p], 

Y t ^ 1 - 5 - m5 2 - aB5 + W t , 

with 

m = sup \b(z)\. (4-31) 

s£ z ^ 1 



Therefore, for m5 < 1, 



Y t ^l-26-aBS + Wt, te[Q,p], 
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so that 



{F p = 1}d| sup W t > (2 + aB)<?Jn j Q< inf ^ W t > (2 + aB- L)Sj. (4.32) 



Choosing L = 3 + a£> and applying a scaling argument, we deduce that 



P 



s: t s: <5 2 



sup W t > (2 + a£)5 [• n <^ inf W t > (2 + aB - L)5 



P 



sg t sg 1 



sup W t > {2 + aB) \ C\{ inf W t > -l\)=: c" E (0, 1). 



sg t sg 1 



(4.33) 



We note that the above quantity c" is independent of 5 and T. Moreover, we deduce 
from fl432] ) that P(F p = 1) ^ c" and therefore, from ( [430| , that 

p(T,y) ^ (1 — c") exp(K5 2 ) sup sup p(£, z), 

zei(L8)te[o,T] 

with X(r) = [1 — r, 1], for r > 0. Choosing 5 small enough such that (1 — 
c")exp(K5 2 ) ^ (1 - c"/2), we obtain 



sup sup p(t,z), y Gl(6). 
zei(LS) te[o,T] 



Modifying c" if necessary (c" being chosen as small as needed), we can summarize 
the above inequality as follows: for 5 ^ c", 



p{T,y) < (1 - c") sup sup p(t,z), y G Z(<S). 
z eX(L<5) te[o,T] 



(4.34) 



We now look at what happens when T ^ 5 2 in (4.30 ). In such a case p in the previous 
argument must be replaced by p AT. Observing that p(T — p AT, Y pA x) = on the 
event {p ^ T} U {F pAT = 1} (since p(0, ■) = on [1 - e/4, 1]) and following fl430] ), 
we obtain, for y e 

p(T,y)< exp(^ 2 )[l-P({F pAT = l}U{p^T})] sup p(M), (4.35) 

(t,*)eS'(«,i) 



with Q'(5,L) = e [0,T] x [1 - Now, th e rig ht-hand side of gjjgp 

is included in the event {I^at = 1} U {p ^ T} so that (4.33) yields a lower bound 
for P({Yp A T = 1} U {p ^ T}). Therefore, we can repeat the previous arguments in 
order to prove that (4.34) also holds when T ^ 5 2 , which means that (4.34) holds 
true in both cases. 

Therefore, by replacing T by t in the left-hand side in (4.34) and by letting t vary 
within [0, T], we have in any case, 

sup sup p(t,y) ^ (1 — c") sup sup p(t,z). 
yeZ{6) te[o,T] zgX(L5) te[o,T] 

By induction, for any integer n ^ 1 such that L"^ ^ r , with r = c" A (e/4), 
sup sup p(t,y) ^ (1 — c") n sup sup p(t,z). 

y£l(S)te[0,T] zel(L"5)te[0,T] 



38 



F. DELARUE, J. INGLIS, S. RUBENTHALER, E. TANRE 



Given 5 G (0,7*o/L), the maximal value for n is n = |_hi[ro/<5]/lnLj . We deduce 
that, for any 5 G (0,r /L), 

(4.36) 



sup sup p(t, y) ^ (1 — c' 
j/£Z(5) te[o,r] 



,//\(lri[ro/5]/lni)-l 



sup sup p(t, z) 

:6Z(e/4) ie[0,T] 



Following (4.16), we know that 



sup sup pit, z) ^ sup sup 

z&l(e/4) te[0,T] zeX(e/4) fe[0,T] 



ct 



cxp 



0? 



£012 



c^t 



(4.37) 



for some constant Ct only depending upon T and K and where (^ )o ^t^T stands 
for the solution of the ODE 

' ' $0 =1 -I'D- 



dt 



b(# t ) +ae'(t), te[0,T\ 



Pay attention that we here use the same notation as in (4.17) for the solution of the 



above ODE but here e(t) is not given as some E(M 4 ). Actually, we feel that there 
is no possible confusion here. Notice also that e is fixed and does not depend upon 
the initial condition xq. 

By the comparison principle for ODEs and then by Gronwall's Lemma, we deduce 
from the fact that e is B- Holder continuous of exponent 1/2 that 

0f <: 0j- e ^ [1 - e + At + Bt 1/2 ] exp(At), t G [0, T\. 



Using the above inequality, we can bound the right-hand side in (4.37). Precisely. 



the above inequality says that the exponential term in the supremum decays expo- 
nentially fast as t tends to so that the term inside the supremum can be bounded 
when t is small; when t is bounded away from 0, the term inside the supremum is 
bounded by ct / y/i. It is plain to deduce that 

(4.38) 



sup sup p(t, z) ^ c-r, 

2£X(e/4) te[0,T] 



for a new value of ct, possibly depending on e as well. Therefore, for 5 G (0,r /L), 

sup sup p(t,y) ^ 



(4.36) yields 



Or 



y&(6) te[0,T] 



with rj = — ln(l — c")/lnL. This proves (4.26) for y G (1 — r^/L,!). Note that r] 



is here independent of T, contrary to what is indicated in the statement of Lemma 



4.7 However, we feel it is simpler to indicate T in r/x as the constant B in the sequel 



will be chosen in terms of T thus making r\ depend on T. Using (4.38), we can easily 
extend the bound to any y G (1 — e/4, 1) by modifying if necessary the parameters 
/it and tit therein. This completes the proof. □ 

4.5. Bound for the gradient. Here is the final step to complete the proof of 
Theorem I4TT1 



Proposition 4.8. For any e G (0,1), any B > and any T > 0, there exists a 
constant M. T > 0, only depending upon T, B, e, K and A, such that, for any initial 
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condition xq < 1 — e, any continuously differentiable non- decreasing deterministic 
mapping [0, T] 3 1 1— > e(t) satisfying 

e(0) = ; e(t) - e(s) < B(t-s) 1/2 , O^s^t^T, 

if (Xt)o < t < t denotes the solution of the SDE 

dXt = b(xt)dt + ae'(t) + dW t , t E [0, T] ; % = x , 

t/ten ; /or any integer n such that [4/e~|, 



1 - expf-Afy (1 + aC T )n 



— (l + aC T ), te[0,T], 



where p(t : y) is the density of Xt killed at 1 as in (3.3), r/y zs as in Lemma J^.l, and 

Ct = sup e'(t). 

o s£ t ^ T 

Proof. We consider the barrier function 

q{t, y) = 6 exp(Kt) [l - exp( 7 (y - 1))] , t^O, yER, (4.39) 
where 7 and G are free nonnegative parameters. Then, for t > and y < 1, 

+ (b(y) + ae'(t))d y q(t,y) - ^d 2 y q(t,y) 
= exp(irt) exp( 7 ( 2/ - 1)) (-(%) + ae'(t)) 7 + ^7 2 ) + Kq(t, y). 



Keeping in mind that 



sup e'(t) = Ct, 

< t < T 



and choosing 

7 = 2(max(m, 1) + aC T ), (4.40) 
where m = sup < z < x as before, we obtain, for t E [0, T] and ?/ E (0, 1), 

-(%) + ae'(t))7+ ]p 2 ^ - 2(max(m,l) + aC r ) 2 + 2(max(m, 1) + aC T ) 2 = 0. 

Thus, for t E [0,T] and y E (0, 1), 

d t q(t,y) + (b(y) + ae'(t))d y q(t,y) - \d 2 q{t,y) > Kq(t,y) ^ - b'(y)q(t,y), 



which reads 

d t q(t,y) + d y [(b(y) +ae\t))q(t,y)] - \d 2 yy q{t,y) ^ 0. 

For a given integer n ^ [4/e] , we choose as the solution of 

2(max(m, 1) + aC?) 







1 — exp 



n 



-'IT 



(4.41) 



(4.42) 
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with [It and rjT as in the statement of Lemma 4/7 Pay attention that the factor in 
the left-hand side cannot be as max(m, 1) > 0. Notice also q thus depends upon 
n. By Lemma 4/7 , we deduce that 

1 



n 



n 



Now, we can apply the comparison principle for PDEs. Indeed, we also observe that 
q(0,y) ^ p(0, y) = for y G [1 - 1/n, 1] and 1) = p(t, 1) = for t G [0,T]. 
Therefore, by (4.41), we have 



p{t,y)<q(t,y), te[0,T), y e 



l 



n 



1 



Since 1) = = q(t, 1), we deduce 

2/i T (max(m, 1) + aCy)n"' ,T 



1 — exp[— 2(max(m, 1) + aCV)/n] 



(4.43) 



exp(ATi). (4.44) 



□ 



We now complete the proof of Theorem 4.1 We make use of Proposition 3.4 



Recall (3.6) 



e\t) = - ^ \d yP ^ s \t - s, l)e'(s)ds - l -d y p(t, 1), t G [0, T], 

where p represents the density of the process X killed at 1 and p(°' s *> represents the 
density of the process X^ s driven by e" s = e(- + s) — e(s) (see ( 4.12[ )) killed at 1 with 



Xq = as initial condition. 



By Proposition 4.5 and Lemma 4.3, we know that, for a given s G [0,T) and 
for the prescribed values of a, the mapping [0, T — s] 3 r h- > e" s (r) is l/2-H6lder 
con tinuous, the Holder constant only depending upon T, a, e, K and A ( Prop osition 
4.5 permits to bound the increments of e" s on small intervals and Lemma 4.3 gives a 



trivial bound for the increments of e" s on large intervals). Therefore, by Proposition 

M T n- r > T 



4.8 we know that 
\d y p^(t-s,l)\^ 



I ir'nx7M^ (1 + a ^' te[s,T], (4.45) 
1 — exp[— M T (1 + aCrjn 1 \ 

for n ^ [4/e] and for some constant M.? only depending upon T, a, e, K and A. 
The same bound also holds true for d y p(t, 1). 

We deduce that, for any t G [0, T] and any n such that n ^ [4/e] , 



e'(t) < 



Vt 



1 - exp[-A4y X (l + aCV)™- 1 ] 



1 + aC T 



e(T) + 1 



By Lemma |4/3 we have a bound for e(T) = E(M^), which means that we can bound 
(e(T) + l)/2 in the right-hand side above by modifying the constant M.?- Recalling 



Ci 



sup e'(t), 

o t sg T 
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we deduce that 



C T (l - exp[-A^ T 1 (l + aC^rT 1 ]) ^ M T (1 + aC T ) 



n 



-VT 



(4.46) 



Choosing n large enough such that the right hand side is less than (1 + o;Ct)/2 (so 
that n depends on T) and multiplying by a, we get (since a G (0, 1)): 



olCt 1 1 
— ^ - + aC T exp[-M T \l +aC T )n 



-ii 



This shows that olCt must be bounded in terms of M.t and n. Precisely, we have 
aC T ^ 1 + 2 sup [rexpl-M? 1 ^ + r)n' 1 }] := R < +oo. 



r > 



By (4.46), we deduce that 



Ct ^ sup 

^ r s£ R 



M T (l + r)n' riT 



1 - exp[-MT 1 (l + r)n- 1 ]_ 
which is independent of a (for a G (0,o;o]) ? as required. 



□ 



5. Existence and uniqueness for all time 

We now put everything together to arrive at our goal, which is Theorem 5.2 below. 

Lemma 5.1. For any T > 0, initial condition Xq — Xq < 1, and a < ciq, where 
cto = 0^0(^0) i> s as i> n Theorem 4-F there exists a constant Cdcn{T) depending only 
on T, xq, K and A such that any solution to (2.1) on [0, T] according to Definition 

d 



2.3 satisfies 



dy 



F(X t edy) ^C dcn (T)(l-y), 



for all y G (1 — e/8, 1) and t G [0, T], with e = min(l, 1 — xq) 



Proof. We assume that (X t )o ^t^T is the unique solution to (2.1 ) with Xq = xq up 
until time T, and set e{t) = E(M 4 ). Following the notation of Section [3] (see also 
the last part of the proof of Theorem 4.1), for y ^ 1 and t ^ T, let 

p(t,y):=j- W{X t edy,t<n), 

P ( °' s) (t,y) := (x}' edy,t< 4°\Xt = 0) . 



By Theorem 4.1, we know that e is A^r-Lipschitz continuous, so that by (|4.43|), 

p{t,y)^q{t,y), te[0,T],ye 



n 



where n stands for [4/e] and q is given by (4.39), with 7 and O being fixed by (4.40) 
and (4.42), with Ct = M.t- By the specific form of g, this says that there exists a 
constant C T , depending only on T, xq, K and A, such that 



p(t,y)^C' T (l-y), te[0,T), ye 



1 - 
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using the elementary inequality 1 — exp( 
argument applies to p(°' s \t — s,y), i.e. 



-x) 



< x for x G 



Clearly, the same 



P 



(0,8) 



(t-s,y)^C' T (l-y), 0O<KT,|/6 



Now, following the proof of (4.11), we get for t G [0,T] and y G [1 — e/8, 1], 

j-F(X t G dy) = pit, y) + [ t p^ s \t-s,y)e'(s)ds 

ay Jo (5.1) 

<C r (l + e(T))(l-y), 



where we use Lemma 3.3 for justifying the passage to the density in (4.11). By 
Lemma 4.3, this completes the proof. □ 



Theorem 5.2. For any T > 0, initial condition Xq = Xq < 1, and a < ao, where 
«o = cto{xo) is as in Theorem there exists an unique solution to the nonlinear 
equation ( 2.1[ ) on [0,T] according to Definition 2.3 



Proof. We would like a solution up until fixed time T > 0. The idea is to iterate 
Theorem 3.1, which is possible thanks to Lemma 5.1 Indeed, by Theorem 3.1, we 



have that there exists a solution to (2.1) with Xq = xq up until some small time 



T\ > 0. By Lemma [5. 1[ we therefore have that 
d 



dy 



F(X Tl edy)^C den (T 1 )(l-y), ye 



1 - 



where e = min(l — Xq, 1). If T\ ^ T we are done. If not, we have the above density 
bound for (d/ dy)F(Xx 1 G dy) and we have an explicit bound for EjX^ | deriving from 



Lemma 4.3 and from Gronwall's Lemma. We also know from (5.1) and Lemma 3.3 
that the density of Xt x is different iable at y — 1. Therefore, we can apply Theorem 



3.1 again to see that there exists a solution to (2.1) on some interval \T\ } T\ + T2] 
starting from X^. As T2 depends upon X^ through e, Cden{Ti) and E|XrJ only 
and as these quantities can be bounded in terms of T, e, K, A only, we then see that 

T 2 ^ 0(T) 

for some constant <f)(T) that refers to T, a, e, K, A only. Now we know that there 
exists a solution to (2.1 ) with X = x 011 [0> T\ + T 2 ]. If Ti + T 2 > T we are done. 



If not, by Lemma [5.1 1 once again 
d 



dy 



P(A Tl+T2 G dy) ^ Cde^Ti + T 2 )(l -y), ye 



1 



and we can then repeat the argument n times to get a solution up until time T\ + 
• • • + T n , where all Tk ^ 4>{T) for k ^ 2 i.e. each time step is of size at least <fi(T). 
It is then clear that there exists n ^ 1 such that T\ + • • • + T n ^ T, and so we are 
done for the existence of a solution. 

Uniqueness of the solution proceeds in the same way. Given another solution 
(A t ', M f ') ^ t ^ t on the interval [0,T] in the sense of Definition 2.3, it must satisfy 



the a priori estimates in the statements of Theorem |4.1| and Lemmas |4.3| and |5.1[ 
In particular, dividing the interval [0, T] into subintervals of length <p(T) (except for 
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the last interval the length of which might be less than <p(T)), with the same <f>(T) 
as above, we can apply the contraction property in Theorem |3.1| on each subinterval 
iteratively. Precisely, choosing A\ accordingly in Theorem 3.1[ we prove by induction 
that the two solutions coincide on [0, 0(T)], [O,20(T)], and so on. □ 

Remark 5.3 (A closed formula for e(t) = E(M t )). This remark is motivated by the 
proof of finite time blow-up in [3] (Theorem 2.2). In fact, thanks to the knowledge 



about the existence of a solution for a sufficiently small provided by Theorem 5.2 , 
we can extract some interesting formulae by modifying the methods in [3]. 

Indeed, fix Xq < 1 and suppose that we are in the Ornstein- Uhlenbeck case, so that 
b(x) = —Xx with A ^ 0. Let a = «o(^o) be as in Theorem 5.2. Then for a < a® a 
solution to (2.1) with X = Xq < 1 exists for all time in the sense of Definition 2.3. 
Set e(t) = E{M t ). 

For jj, > 0, let (p(x) = exp(fix). Then since e(t) is continuously differ entiable, by 
the generalized ltd formula we see that for all t ^ 0, 

E[<p(exp(Xt)X t )\ = E(ip(X )) + f aj2exp(Xs)e'(s)E[ip(exp(Xs)X s )]ds 



+ 



+ 



ir 



exp(2Xs)E[(p(exp(Xs)X s )]ds 
[<p(0) — ip(exp(Xs))]e'(s)ds. 



Differentiating with respect to t and writing z(t) = E[<p(exp(Xt)X t )] yields 



d 
dt 



z(t) 



afiexp(Xt)e'(t) + — exp(2At) 



z(t) + [<p(0) - p(exp(At))]e'(*) 



for allt > 0, with initial condition z(0) = exp(fiXQ). Writing r](t) := afiexp(Xt)e'(t)-\- 
(/x 2 /2) exp(2At) andv{t) := [ip(exp(Xt)) — Lp(0)]e'(t) ^ ; we can solve this first-order 
ordinary differential equation in the standard way to see that 



z(t) = exp ( / rj(s)ds 



v(s) exp ( — / r/(u)du I ds 



3(0) - 

Now, by the definition of ip, we can note that z(0) ^ ; and that 

z(t) exp(— /x exp (At)) ^1, t ^ 0, 

since X t < 1. Setting 



7(f) := exp I — /i exp (At) 



r)(s)ds 



it follows that 



^ 7 (t) 



<0)- 



v{s) exp 



r)(u)du ) ds 



^ 1, 



(5.2) 
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for all t ^ 0. But we can see that ^(t) — > oo as t — )■ oo /or all fi > 0. Indeed, 
7(t) = exp ^— /iexp(At) + ^ [a//exp(As)e'(s) + (/i 2 /2) exp(2As)] ofs 
^ exp f — /iexp(At) + — J exp(2Xs)ds 



exp ( /iexp(At) -^-exp(At) — 1 
L4A 



exp ( -/i + y t 



EL 
4A 




oo. 



Thus for the bound (5.2) to remain true when t — > oo, we must have 

z(0) — J v(s) exp ^— J r](u)du^j ds = 0, V/i > 0, 

which provides a surprising closed formula for e' . We can then perform an integra- 
tion by parts. Indeed, by writing e(s) = J Q S exp(Xu)e'(u)du for all s ^ 0, the above 
equality yields 

f°° ( \ 

= z(0) — J [</?(exp(As)) — tp(0)]e'(s) exp I — afie(s) — — [exp(2As) — 1] J ds 



z(0) 



0(s)afj 1 e' (s) exp (— afie(s)) ds, 



with the convention that [exp(2As) — 1]/(4A) = s/2 if X = 0, and where 



^(exp(As)) - <p(0) 



Since it is clear that 0(oo) = 0, we can perform an integration by parts, so that 

POO 

= z(0) — 0(0) — / v'(s) exp (— afie(s)) ds. 
Jo 

In other words, whenever a < ato, we must have that the map t i-> e(t) = E,(M t ) 
satisfies 

ef - 1 



exp(/iXo) 



OL[L 



+ / v'(s) exp (— a\ie(s)) ds, 



(5.3) 



/or all fi > 0. It seems to be an interesting question to determine whether this 
characterizes e and more generally to determine whether the existence of a solution 
to (5.3) implies the existence of a global solution to (2.1). In particular, in the case 

(5.4) 



when X = 0, (5.3) reads 



E 

2 



f ( , x A* 2 A a//exp(//x ) 

/ exp -afxe(s) - — s ) ds = 1 — — . 

Jo \ 2/ exp - 1 
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6. Explicit bounds for the criticality 



Theorem 5.2 says that, for a given initial condition Xq, existence and uniqueness 



of a solution to ( 2.1[ ) holds when the self-excitation rate a lies in (0, ao(xo)), for 



some ao(xo) G (0, 1], the value of which is given in the statements of both Theorem 



4.1| and Proposition 4.5 This raises several questions: 

(1) Does the model indeed exhibit several regimes for global solvability according 
to the values of of! 

(2) If so, are there two phases exactly (using an analogy with statistical mechan- 
ics)? To put it differently, does there exist a critical threshold a c such that 
global solvability holds for a < a c and fails for a > a c l Or are there any 
more complex cases when global solvability fails for some a + and then holds 
for some a~ > a + l 

(3) Conversely, are there any cases when global solvability holds for any a < 1? 



Although Theorem 5.2 does not provide any answers to these questions, it does say 



that, for a given initial condition Xq < 1, the critical rate 

a~(x ) = sup{a _ G (0,1] : Q(a,x ) holds for any a G (0,a~)} 
is (strictly) positive, since it is larger than a (x ), where Q(a, x ) stands for the prop- 



erty of global unique solvability of equation (2.1) when driven by a and initialized 
with Xq. 

In this framework, the lower bound a~(xo) ^ ao(xo), with ao(xo) given by Propo- 
sition 4.5 and Theorem 4.1, is complementary to the previous result in |3]: there it 
is indeed proven that the existence of a global solution fails for some values of a < 1 
and some initial conditions, including Dirac masses located right below the barrier 
1, in the case when b(x) + Xx ^ 0, with A > 0. This shows that, in some cases, 
a~(xo) is indeed strictly less than 1. Actually, the result in [3] says more: for any 
value of A > 0, there exists a critical threshold x c < 1 such that, for any Xq G (x c , 1), 
there exists some a' (xo) < 1 such that global existence fails (in the sense that e' 
blows up in finite time) for any a G (a' (xo), 1). In other words, the critical rate 

a^(xo) = inf{a + G (0, 1] : Q(a,Xo) fails for any a G (ct + , 1)} 

is strictly less than 1 for xq close to 1, with the convention that a^(xo) = 1 if the 
set in the right-hand side is empty. Obviously, 

a~(x ) ^ a+(x ), 

but, as written above, the following questions remain open: 

(1) Are there any cases when a~(xo) ^ o^{xq)1 

(2) Are there any cases when a~(xo) = 1? 

More generally, finding the optimal values of a~(xo) or/and a^(xo) turns out to be 
a really challenging problem. 

In this section, we address two problems: (i) the first is to provide an explicit 
lower bound for a~(0) in the case when b(x) = —Xx, x G M, for some A > 0; (ii) 
the second is to provide a lower bound for a~(xo) and an upper bound for a^(xo) 
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in the case when 6 = but xq is chosen in such a way that a+(xo) is known to be 
less than 1 and then to discuss the bounds in light of some numerical experiments. 

6.1. Ornstein-Uhlenbeck case. In this paragraph, we investigate further the case 
when b(x) = —Ax (referred to as the Ornstein-Uhlenbeck case, since the dynamics 
of X up until the first jump time coincide with those of the Ornstein-Uhlenbeck 
process), which is the choice most often made in neuroscience. As above, we assume 
that the diffusion coefficient is 1. We then exhibit a lower bound for the optimal 
value of a~ according to the methodology described in Section |4j As the method 
involves rather tedious computations, we will focus on the case X Q = only: this 
makes things a bit simpler. We thus consider 



X t — —A / X s ds + ae{t) + W t -M t , 0, 
Jo 



(6.1) 



where e(t) = E(Mf). When Xq = 0, as is the case here, numerical experiments 
indicate that existence and uniqueness of a global solution may hold for any a G 



(0, 1). Indeed, equation (6.1) is then solved by using a particle method: we refer to 



the forthcoming paper [5] for a complete description and analysis of the numerical 
scheme. For instance, when A = 0.01, we get as a plot of the mapping t (->■ e(t) 
associated to the approximate solution that given in Figure [TJ 




Figure 1. Plot of t h> e{t) for X = 0, A = 0.01, a = 0.9 and t £ [0, 10]. 



Numerically speaking, existence and uniqueness is then suspected to hold since 
the graph of the mapping t h- > e(t) is smooth and, in particular, has no jump. We 
refer to Subsection 62 for a counter-example when a jump occurs for some xq ^ 0, 
indicating that existence of a global solution, as defined above, cannot hold. Below, 
we are not able to prove that a~(0) = 1, showing that our methodology is not 
optimal. Anyhow, we manage to give an explicit lower bound for the critical value 
of a~(0) and to prove that it converges toward 1 when A tends to +oo. 

The strategy is to seek ao such that, for any t ^ and y in a neighbourhood of 



1, [d/dy]P(X t G dy) < 1/a for all a < ao, so that the condition (4.5) in Lemma 4.4 
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is satisfied (which then allows us to pass to an unique global solution thanks to the 
results of Sections 111 and IBT). 



First Step. The first step is to provide some a priori estimates for the solution 
on any interval [0,T] where existence holds. Precisely, the point is to refine the a 
priori L°° bound for the function (e(i) = E(M t ))o <g t < r given in the statement of 
Lemma 4^3, by taking advantage of the mean-reverting property of the Ornstein- 
Uhlenbeck process. As in Section |4| in all the a priori estimates below we assume 



that (X t , M t )o ^t^r stands for a solution to (6.1 ) in the sense of Definition 2.3 
Here is the first claim: 

Lemma 6.1. For any t G [0,T], 

:(E[[(X,)-] 2 ] 1/2 , E[(X t )+\, (1 - a)exp(-Xt) j\xp(Xs)e'(s)ds^ </ a (A,f), 

with fi(X,t) 



max 



1 - exp(-2At) 



2A 



where (a;)_ = max(— x, 0) for x G 



1/2 



Proof. We first investigate how far the process (X t ) Q ^ t can go in the negative 
direction. The point is to apply Ito's formula to ([(X t )_] 2 )o ^ t ^ t' 

[(Xt)-} 2 

= 2A / (X s _)_X s ds-2a [ (X 8 _)_e'(s)ds + [ l {Xs _ <0} ds-2 [ (X s _)_dW s 
Jo Jo Jo Jo 

= -2A / (X s ) 2 _ds + 2a [ X s l {Xs<0} e'(s)ds+ [ l {Xs<oy ds-2 [ (X s _)_dW S} 
Jo Jo Jo Jo 

for t G [0, T\. Here, stands for the negative part of the left-limit of X at time 

s. To pass from the second to the third line, notice that X s „ = X s when X s _ ^ 0. 
Multiplying by exp(2At), reapplying Ito's formula and then taking the expectation, 
we deduce: 



E[exp(2At)(X t ) 2 _] 



o 



* exp(2At) 

exp(2As)as = — — 

2A 



exp(2At)/ 2 (A,t). (6.2) 



Now, multiplying in the same way (6.1) by exp(At) and applying Ito's formula, we 
get: 

exp(Xt)X t = a exp(Xs)e'(s)ds + / exp(Xs)dW s — / exp(As)<iM s . 
Jo Jo Jo 

Taking the expectation, we deduce: 

I - n ) / exp(Xs)e'(s)ds = -E[exp(\t)X t ], 
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so that, writing X t = (Xt)+ — (X t )-, with (x 
E[exp(At)(X t )+] +(!-« 



= max(x, 0), 

exp(As)e'(s)ds = E[exp(At)(X t )_] 
^ exp(A£)A(A,t), 



where the last inequality follows from (6.2) and Holder's inequality. 
We deduce: 



(6.3) 
□ 



Lemma 6.2. For any t,t' G [0,T], t < t' 



e(t') - e(t) ^ 



f 2 (X,t'-t,T) 



a 



with 



/ 2 (A, r, T) = min [exp^^A, T) (l + Ar) , 1 + (A/2) 1 ^ + (2r/vr) 1 /2] . 

Proof. We prove first that the bound holds with respect to the first argument in the 
minimum defining f 2 . By Lemma 6A_ , we know that, for ^ t ^ t + h ^ T, 

rt+h fi{X,T) 



r f ( 

[e(t + h)- e(t)) exp(-Xh) ^ exp[-A(t + h)) / exp(Xs)e'(s)ds ^ — 

Jo 1 

Choosing h = 1/A, we deduce that for < t < s ^ T A (t + 1/A) 

,A(A,T) 



a: 



e(s) - e(i) ^ exp(l)- 



Summing over a subdivision of [t, t'} made of small intervals of lengths less than 1/A, 
we complete the proof for the first argument in the minimum. 

We now turn to the second argument in the minimum. As in the proof of Lemma 



4.4[ we can assume without any loss of generality that t = by shifting the system 
(the shift being of length t), provided the new initial condition X is chosen as an 
arbitrary random variable. Precisely, the new initial condition X Q is then understood 
as the previous X t . Similarly, the new variables (X s )o t '-t are understood as the 
previous (X s ) t <: s <: t' so that, by Lemma 6.1 , we can assume the new (X s )o <g s <g t >-t 
satisfy 



e [(*,)+] </!(A,r)< 



2A 



1/2 



Following the proof of (2.10)— (2.11 ), we then get for any s G [0,f — t\. 



e(s) = E(M S ) < E 



sup (Z r ). 

< r < s 



< 1 + A 

< 1 + 



E[(X s ) + ]ds + ae(s) +E 

1/2 



sup W r 

< r < s 



(f -t) + ae(s) + E 



sup W r 

< r < s 



so that, for any s G [0, if — t], 



els) ^ 



1 + (A/2) 1 / 2 ^ -t) + y/2(t'-t)/n 



1 — a 
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the last term in the right-hand side following from 



E ( sup W s 

K < s < t 



(6.4) 



Indeed, by the reflection principle, 



P ( sup W s ^ a ) = 2P (W t ^ a) - 1 

< t 



- z2 ' 2t dz, 



so that 



E ( sup W s 

kS < t 



ze 



- z2 ' 2t dz 



Second Step. According to the program we fixed in the proof of Lemma 4.6 



□ 



we 



now investigate the bound of the density in small time. In the Ornstein-Uhlenbeck 

m 

1 



case, (4.16) (with x = 0) has the form: 
d 



-P ( Z? e dx 
dx ' ' 



x exp 



(2 7 r)V2/ 1 (A,t) 



2A 2 (A,t) 



x — aexp(— At) / exp(Xr)e'(r)dr 



(6.5) 



for t ^ T, with f\(X,t) as in the statement of Lemma 6.1 and Z t ° given by (4.14) . 
More generally, ( |4.18[ ) reads: 



P ( zh° e dx 



d 

dx 
x exp 



(2 7 r)V2/ 1 (A,t- s ) 



2f!(X,t-s) 



t-s 



x — a exp[— X(t — s)] / exp(Ar)e'(s + r)dr 



2-1 



for s < t ^ T. By Lemma |6.1[ we notice that the "mean-term" in the exponential is 
bounded by 



a i 



t-s 



exp[— A(t — s)] / exp(Ar)e'(s + r)dr = a exp(— At) / exp(Ar)e'(r)(ir 



a/iCM) 
1 — a 



(6.6) 



By Lemma 6J2 it is also bounded by 
exp [— A(£ — s 



a 



3 exp(Ar)e'(s + r)dr < a(e(t) - e(s)) s= (6.7) 

1 — a 



Therefore, since r i-> f2(X,r,t) is non-decreasing, if, for some e G (0, 1), 



a 

1 — a 



mm 



[/i(A,t),/ 2 (A,t,t)] ^1-e, 
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then, for x G (1 — e5, 1), with 5 G (0, 1), we have for any s G [0,t) (pay attention 
that this covers ( 6.5[ ) as well): 



d 

dx 



P ( G cfe ) ^ 



(2 7 r) 1 /2 /l (A ) t- 

'l-8)e 
(1 - V 



exp 



1 - 5) 2 e 2 



where 



9{u) 



sup 

v>u 



(2tt) 1 /2 



exp 



2\ n 



(2tt) 1 /2 

1 



exp 
exp 



~2 
1 
2 



if m ^ 1, 



if it < 1, 



[ (2tt)V2 

using the fact that u i— >• g(u) is non-increasing and r i— )■ fi(X,r) is non-decreasing. 
Therefore, by ( |4.15[ ), for x G (1 - 1) 

Define T = 1 A A -1 . Then, by Lemma 6.2, for t ^ T A T , 



d_ 
dx 



F(X t edx) < 1 + 



/ 2 (A,T ,T ) 



1_ Q ; (l-5)e* V/i(A,T ) 



again since g is non-increasing, provided 



a: 



min[/ 1 (A,T ),/ 2 (A,T ,T )] <l-e 



(6.8) 
(6.9) 



holds true. 

Third Step. We now investigate the bound of the density for t ^ To. We first 
discuss tti defined by (4.23). By (4.25) (plugging inside the current value of To), and 
the same arguments as in the second step, 



7Ti < 



1 



-9 



[l-6)e 



(e(t)-e(t-T )), 



provided 



a 



1 — a 



l-5)e" VA(A,T ). 
min[/i(A,oo),/ 2 (A,T ,oo)] ^ 1 - e. 



(6.10) 



Pay attention that, in /i(A, •) and /2(A,To, ■), • is understood as To in (6.9) whereas 
it is understood as oo in (6.10). The difference between them comes from the fact 
ly be 



that in this step t may be large. By Lemma |6.2 

1 



l-S)e 9 \f 1 (\,T ) 



l-5)e\ / 2 (A,T ,oo) 



1 — a 



(6.11) 



It now remains to investigate 7r 2 also defined in (4.23). We start with the case A ^ 1, 
so that To = A -1 . Following (4.24) and (6.5) (keeping in mind that To A = 1), we 
have 

7T 2 < SUp 7T 2 (z), (6-12) 
z s£ 1 
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with 

7r 2 (» 



(2tt)V2 /i (A,1/A) 

1 



x exp 



x — exp(— 1) 



z + a exp(Xr)e'(t — T + r)dr 
Jo 



2/! 2 (A,l/A) 

Actually, the bound can be refined for A large by taking into account the fact that z 
is understood as a realization of X t -T and that the probability that X t -T is away 
from is very small when A is large because of the mean-reverting property of the 
Ornstein-Uhlenbeck process. Precisely, we claim: 



7r 2 < max ( sup tt 2 (z), sup -n 2 {z)¥(X t ^ TQ ^ 1/2) 

k z < 1/2 z < 1 

^ max I sup 7r 2 (-2), (2/A) ^ 2 sup 7r 2 (,2) 

z< 1/2 z < 1 



(6.13) 



the last line following from Lemma 6.1 and Markov's inequality. As above (see (6.6) 



and (6.7)), we know that the mean term in the exponential in 7r 2 (z) is bounded by 

exp(— l)z + aexp(— 1) / exp(Ar)e'(t — T + r)dr 

Jo 

fib 

= exp(— l)z + aexp(— AT ) / exp(Ar)e'(t — T + r)dr 

Jo 

= exp(— l)z + aexp(— At) / exp(Ar)e'(r)dr 

Jt-Tb 



^ exp(— i)z + 



a 



1 — a 



min [f^X, oo), / 2 (A, T , oo)] . 



Therefore, if ( 6.10[ ) holds for e G (exp(— 1), 1), then, for x E (1 — (e — exp(— 1))S, I) 
we get from (6.12) and (6.13): 



with 

f4(S, A) 

f 5 (5, A) = max 
Finally, we have 



1 



tt 2 < min [/ 4 (5, A), f 5 (6, A)] , 

(l-<5) 2 (e-exp(-l)) 2 



(27r)V2 /l( A,l/A) 
1 



exp 



(27r)V2 /l( A,l/A) 



exp 



2A 2 (A, 1/A) 
( l_ 5) 2 (£ _ exp( _ 1)/2 )2n 2 1 / 2 / 4 ((5 l A) 



2/ 1 2 (A,l/A) 



A 1 / 2 



1 — a 



(6.14) 



+ min[/ 4 (5,A),/ 5 (5,A)]. 



provided (6.10) holds true, since [d/da;]P(JQ G da;) < 7Ti + 7r 2 by (4.24) 
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We now turn to the case A < 1. Then, Tq = 1 and 



vr 2 ^ 



1 



(2tt)V2 /i (A,1)- 



Therefore 
d 



9 



1-S)e\ A(A,T ,oc) 



/i(A,T 



+ 



a 



(2 7 r)V2/ 1 (A,i; 



(6.15) 



Conclusion. For our purpose, it is worth noting that S can be chosen as small as 
desired, so that the critical value a~ (0) can be obtained by putting 5 = in all the 



previous work. The reason is quite simple: the strategy for proving Proposition 4.5 



applies provided the estimate of the density (4.5) holds true on a small interval of 
the form (1 — 5, 1), 5 being possibly really small. Obviously, when S tends to 0, the 
resulting Holder constant in Proposition 4J3 blows up, but this does not affect the 
determination of ft~(0). 

Start with the case A ^ 1. From (6.8), ( 6.10[ ) (which implies (6.9)) and (6.14), we 
get as critical conditions for a in order to satisfy (4.5): 



1 + 



min [/i(A, oo), A (A, T , oo 
1 



< 1, 



A(A,T 



a 



f 2 (X } T 0} T )g 



A (A, T 



< 1, 



(6.16) 



a 



min[/ 4 (<y, A), f 5 (6, X)]a+ -/ 2 (A, T , 00)5- 



A (A, T 



a 



< 1. 



a 



When A < 1, we get similar conditions as in (6.16), but the third line in (6.16) is 
replaced by 

1 1 



(2 7 r)V 2/l (A,T ) a + F /2(A ' T0 '° O) ^ 



A (A, T 



1 — ft 



< 1. 



(6.17) 



The first inequality in (6.16) is explicit. The two last (and (6.17) as well) are of the 
form 

ft 2 

cift + c 2 - < 1 <^ c%a — (1 + c% + c 2 )ft + 1 > 



' 1 — ft 



ft < 



1 + Ci + Ca - [(1 + Ci + c 2 ) 2 - A Cl } 1 ' 2 
2ci 



where we used ft G (0, 1) in the two equivalences. Therefore, a lower bound for 
the critical value of a can be computed numerically: we denote it by 01$ as the 
coefficients in ( 6.16[ ) depend upon e G (exp(— 1), 1). Optimizing over e, we get a 
new bound, which must be read as the ft (or fto(0) if indicating the initial condition) 
in Proposition 4.5 and in Theorem 4.1 See Figure [2] for a plot of fto(0) in terms of A 
(the optimization of a £ Q with respect to e being performed numerically over a mesh 
of step size 0.01 for the parameter e) and Figure [3] for a zoom on the initial values. 

The minimum is a nin (0) ~ 0.301. It is reached at A ~ 0.23, which is the point at 
which Figure [3] has a singular behaviour. Actually, such a singularity is expected to 
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A 



Figure 2. Plot of a (0) in terms of A e [0; 80]. 




FIGURE 3. Plot of a (0) in terms of A G [0; 10]. 



be of computational essence and certainly indicates that the computations we have 
performed are far from being optimal, especially for A small. Anyhow, for A large, 
Figure [2] shows that the critical value ct~(0) is (at worse) close to 1, which sounds 
as an interesting result from a practical point of view. Indeed, by scaling, the case 
A 3> 1 and a = 1 (with o as in (1.2)) is equivalent to the case A = 1 and ir<l: for 
a macroscopic value of A, the system is thus uniquely solvable for most of the values 
of a in (0,1) when the noise is small. Actually, the convergence of «o(0) toward 1 
as A — > 00 can be proven rigorously, as stated in Proposition |6.3| below. 



Proposition 6.3. In the case Xq = and b(x) = —Xx (with a — 1), the critical 
value cto in Theorem and Proposition 4-5\ can be chosen such that, for A large, 

-1 



a 



1 + 



1 



exp(-l))(2A) 1 /2 



1 



(1 - exp( 
1.12A- 1 / 2 . 



■1))(2A)V2" 



In particular, «o — > 1 as A — > +00. 



Proof. The proof is quite straightforward and follows from the fact that, in (6.16), 



when e is given, the coefficients of a and a /(I — a) in the second and third lines 
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tend to as A tends to the infinity, so that the conditions are automatically satisfied 
for A large. In the first line, the condition is of the form 



a 



1 + 



(l-e)(2A)Va 
Letting e tend to the minimal authorized value exp 



< 1. 



•1), we complete the proof. □ 



The idea of the proof of Proposition 6.3 is thus very simple: the optimal strategy 
when A is large consists in choosing e close to the minimal authorized value exp(— 1). 
Such a strategy is well-checked from a numerical point of view. Indeed, Figure [4] 
provides a plot of the optimal choice for e in terms of A, when optimizing a e Q with 
respect to e (e being chosen amongst a mesh of step size 0.01): 




Figure 4. Plot of the optimal value of s in terms of A G [0; 40]. 



6.2. A case with blow-up. We now investigate the Brownian case (so that 6 = 0). 
We have already proven that, for any initial condition Xq — Xq < 1, a~(xo) > 0. 
On the other hand, by adapting the strategy in [3], we can prove that o^{xq) < 1 
for some x (note that the result in |3J is written for a non-zero A). 

Indeed, from Remark 5.3, we know that, for any solution (X t , M t ) t ^ to (2.1) 
with 6 = 0, (5.4) must be satisfied. If 



exp(/ix ) ^ 



exp(/x) 



(6.18) 



then (5.4) implies that the integral J °° exp(— afie(s)— fi 2 s/2)ds must be non-positive, 
which is clearly a contradiction. In such cases, existence of a solution must fail. 

Minimizing over /j, in (6.18), we deduce that existence fails for a > inf M>0 [[exp(//) — 
l]/[/iexp(/zx )]], so that a^(x ) ^ inf M>0 [[exp(/x) — l]/[/iexp(/ix )]] (if less than 1, 
which might not be the case). As in the previous subsection, this value must be 
compared with numerical experiments. Indeed, by investigating numerically the 
graphs of e(t) = E(M t ) for different values of a, we can detect the emergence of 
some discontinuity and thus give a bound for at£(xo). 
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Choosing x = 0.8, we find inf M> o[[exp(/i) — l]/[/xexp(0.8/x)]] ~ 0.539, so that 
a+(0.8) ^ 0.54. In Figure [5j we numerically observe that the graph of e is regular 
for a = 0.38 but has a jump for a = 0.39. We also observe that the associated 
particle system keeps stable for a = 0.38 but blows up for a = 0.39. From the 
observations we have for other values of a, it seems that, in this framework, global 
solvability fails for a ^ 0.39 and holds for a ^ 0.38. We thus expect that, in this 
framework, a+(0.8) and a~(0.8) coincide, with a common value between 0.38 and 
0.39. Figure [6] below summarizes these facts (when b(x) = and xq = 0.8). For a 
in region D, by the above argument, we know that global uniqueness fails. For a in 
region C, by numerical experiments it seems that global uniqueness also fails, but we 
are unable to prove it. Similarly, we have proven in this article that for a in region 
A, for some ao(0.8) > 0, an unique solution for all time exists, while numerically 
it seems that global solutions also exist for a in region B. We thus suspect that 
a~(0.8) = a+(0.8) lies in the region in between B and C. 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



Figure 5. Plot of t h-> e(t) for x = 0.8, b(x) = 0, a = 0.38 (red) and a = 0.39 (green). 
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FIGURE 6. Critical regions of a G (0, 1), for x = 0.8 and b(x) = 0. 



To complete Figure M we now explicitly find «o(0.8) in the case 6 = under 



consideration, which can be used in Theorem 4.1 and Proposition 4.5 to prove global 
existence and uniqueness according to Theorem 5.2 Precisely, we prove that ao(0.8) 
can be taken as 0.104. 
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In fact we carry out the calculations for a general Xq G (0, 1). As in the previous 
subsection, it suffices to find an ao(xo) > such that in a neighbourhood of 1 



-T-F(X t G dx) < —. 
ax a 



(6.19) 



for all t G [0, T], whenever (X t )o <g t < t is a solution to (2.1) on some [0,T], T > 0, 
with 6 = and a G (0, ao(xo)), i.e. 

X t = x + ae{t) + W t -M t , e(t) = E(M t ), t e [0,T\. 

Suppose that a ^ 1 — rj for rj G (0, 1) with 77 to be chosen later and xq < 1 — a. 
Then consider x G [xq/(1 — a) +e, 1] with e = (1/2) [1 — xq/(1 — a)]. The following 
simple Gaussian estimate will be our starting point: 

P(X t G dx) = — ¥(Z t - M t edx) = J2 l- p ( z t -kedx,M t = k) 



dx 



k=0 



== ^exp 



(2 + — xo — ae(t)) 1 
2t 



(6.20) 



where as usual Z t = X t + M t , (Z t = xq + ae{t) + W t )o ^t^T thus being a drifted 
Brownian motion. 

We will need an a priori estimate of e(t), similar to the one given in Proposition 



2.2 However, in this simple case when b = we can do better. Indeed, using the 



fact that M t = [sup s ^ t (Z s ) + \ (see ( |43| )) and following (j2TT0J) (|2TTT]) , we see that 
e(t) = E(M t ) = E ( [sup(Z s ) + \ 



s < t 



^ xq + ae(t) + E ( sup W s ) ^ x + ae(t) + 



It 

TV 



1/2 



the last inequality following from ( |6.4 ). This yields 

1 



e(*K 



1 — a 



(x + 



(6.21) 



(6.22) 



where C = a/ 2/7r. As in the previous section we look at what happens first in small 
time t. 



Small t: By (6.22), we have that, for any x G [xq/(1 — a) + £, 1], 



x — x — ae(t) ^ x 



Ca 



1 — a 1 — a 



Ca 
1 — a 



(6.23) 



when y/t < e(l - a)/(2Ca). Let i = [er//(2C(l - r/))] 2 . Then since a ^ 1 - 77 we 
have that 



t C t 







2C(r^^ 



e(l — a) 
2Ca 
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so that (6.23) holds for all t ^ to- At this point we define /3 := C(l — r])/r], so that 
V^o — £/(2/3)- From (6.23), we thus have that, for all t G [0,t ], 



exp 



(x — xq — ae(t))~ \ , / 



2t 



2\ / r-2 



and 



(x + k — x Q — ae(t)) 2 \ ( k 2 

VX]){ 2t ^ 6XP 1-2* 



for k > 1. Hence by (|6_20j) 

— P(Xt G dx) ^ -^L= exp 
ax V27rt 



£ 2\ 1 ^ /A; 2 



t G [0,t c 



Note that 



1 

2" 



Moreover 



sup 

t ^ o 



exp , 



-2\ 1 



2 f 1\ 1 

exp 

7T 



27 e' 



and the supremum is attained when i = e 2 /A. 
We thus arrive at 



— ¥(X t G dx) < \ -exp 

dx V 7T 



l\ 1 1 

2h + 2' 



te[o,t ). 



Large t (i.e. t ^ to): For fixed x G [xo/(l — «) + £, 1] and t > define 

iV(t, x) = ( [ae(t) + x Q — x\ ) . 



(6.24) 



Then, by (6.22) again, we have that 



N(t,x) < (ae(t) +x -x)^ < ( ^-Vi + -^x - x J 

/+ \ 1 — a 1 — a / 



^\s + - Vi) <- v 7 * 



1 — a 



1 -a 



(6.25) 



for all £ > 0. Moreover, for t ^ to, 
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1 oo 

=E 



'2-Kt 



exp 



fc=0 

1 



(x + k — x — ae{t)Y 



'2-Kt 



E 



2t 



exp 



k ^ 

fc:x+fc— xq— ae(t) ^ 1 



(x + k — xq — ae{t)f 
2t 



+ 



1 oo 



exp 



fc ^ i 



2t 



^ttjA; GNU{fl}:^l + io + ote(t) - x} + \ 
2irt 2 

^[iV(M) + 2] + i 



+ 



Ca 



2 V 7T £T 



+ 



(6.26) 



using the symbol jj for the cardinal in the third line, and the fact that y/to = ej (2/3) 
and /? = C(l-r])/r)^Ca/(l-a). 

We now consider two cases. 
Case 1: when the RHS of (6.24) is bigger than the RHS of (6.26). This occurs if 

1\ 1 



vr e J2-K 



exp 



^ ^(4 + e)/3 ^ exp 

C — 



<^ 1 - r/ ^ 



We therefore fix rj such that 1 — r\ 



2 + Cexp(l/2)(4 + e)' 
2/(2 + Cexp(l/2)(4 + e)). Thus we have a set 



of constraints on a in order to achieve (6.19), which are (recall (6.24)) 

-i 



a < 




exp( 




a ^ 1 — rj 



2 + Cexp(l/2)(4 + e)' 
subject to a < 1 — xq. Substituting in the value of e, the constraints become 



a < 



a < 



2Cexp(-l/2)(l -a) 1 

1 — a — Xq 2 



/fV) 



2 + (4 + [1 - a - z ]/ [2(1 - a)]) Cexp(l/2) 



(6.27) 
(6.28) 



still subject to a < 1 — xo, with C = \j2pK. The first thing to note is that fi°(a) 
decreases to on the interval a G [0, 1 — xq). By plotting the curves of and 
/f°(a:) ( a £ [0) 1 ~ x o)) f° r different values of Xq, we see that (6.27) is more stringent 
than (6.28) for large values of xq, while the reverse holds for small values of xq. 
Indeed the critical value of xq is approximately 0.51. 
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Hence for xq ^ 0.51, if (6.27) is satisfied subject to a < 1 — Xq, then we achieve 

xo,5o(xq)] where 5o(^o) is 



(6.19). Thus for xq ^ 0.51, we choose «o(^o 
such that 



min 1 



1 = a (x Q ) 



2Cexp(-l/2)(l-So(x )) 



a {x ) 



X 



which solves as a second-order equation. For instance, this yields a (0-8) 
which completes Figure 6 above. 



0.104, 



Now, if xq J$ 0.51 then (6.28) is the more stringent. In this case we must solve 
5o(£o) = /^(So^o)) and then choose ao(xo) = min[l — Xo,5o(a;o)]. For example, 
c^o^o) ~ 0.264 when xq = 0.4. We note here, however, that as xq \. 0, /2°( a ) ~~ 
4/ [4 + 9Cexp(l/2)] « 0.2525 for all a G (0, 1). We can in fact hope to do better 
when xq is small. 

Case 2: when the RHS of (6.26) is bigger than the RHS of (6.24). In order to 



try and do better for x small, we now look at this case, which occurs when (4 + 
e)B ^ 2exp(— 1/2): here e depends upon Xq but B does not, so that (B acts as a free 
parameter that can be tuned accordingly to the value of Xq. Such a strategy turns 
out to be efficient when xq is small. Indeed, our constraints now become 



a < 



a < 1 




B+C B+C 



(6.29) 
(6.30) 



again subject to a < 1— xq, where we may vary /3, as long as (A + e)(3 ^ 2exp(— 1/2). 



Now, optimizing (6.29) and (6.30) over B by equalizing the two right-hand sides and 
thus by solving 



+ 1 



CB 1 

— + - 
2 2 



(B + C & 



4 



1 \C(B 2 - B - 2C = 0, 



we would like to take B such that 



2( - + 1 )CB = l + x \l 



I | ^ + 1 ) 2C.'-'. 



(6.31) 




This can be done provided it remains true that (4 + e)B ^ 2exp(— 1/2), i.e. that 

' 4 + i' 

which we can see numerically holds for e > 0.15 (which basically says that xq must 
be small, as announced). 

If we suppose that a ^ 0.42 (we will see below that we do not in fact lose 
anything by doing this), then since e — (1 — xq — ct)/2(l — a), it is certainly true 
that e > 0.15 when xq is small. Indeed it is true when xq < 0.28. In this case, 
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a ^ 0.42 < 1 - 0.28 < 1 

a (x ) = 



Xq. By (6.29) and (6.31), we then take «o(^o) such that 
-l 




-i 



1 - a (x ) - x 

Then it follows that whenever xq < 0.28, ( |6.19[ ) holds for all a < ao(^o)- Taking 



xq = we calculate that 

«o(0) 



4 ( 3 + v 7 ! + 72C 2 \ 



0.41, 



which is much better than that achieved above for small Xq. This must be also 
compared with the value obtained for «o(0) in the Ornstein-Uhlenbeck setting when 
A = 0: referring to Figure [3j it is approximately equal to 0.37, which is a bit less 
but not so different. The difference between both is not so surprising. Indeed the 
method used in the Brownian setting is more straightforward since it benefits from 
the Gaussian property of Z t . 



7. Appendix: Proof of Lemma l3~3l 



First Step. We first discuss the solvability of (3.4). We start with the following 
case: we assume that xo — x o £ ( — °°> 1)> an d that b and e are smooth and bounded, 
with bounded derivatives of any order. Then, by |U Th 1.10, Chap. VI], we know 
that the generator of Xi namely the family of second-order differential operators 

(£*,,:= {b{-)+ae'{s))d x + \dl x ) Q<s<T , 

admits a Green function G : [0,T] 2 x (— oo, l] 2 3 (s,t,x,y) h-> G(s,x,t,y). For a 
given (t, y) G [0, T] x (— oo, 1], the function [0,t) x (— oo, 1] 3 (s,x) h> G(s,x,t,y) 
is a classical solution of the PDE 

d s G(s, x, t, y) + C S)X G(s, x, t, y) = 0, 

with G(s, l,t, y) = 0, for s G [0,t) and G(s,x,t,y) — > 5 (x — y) as s /* t, where 
Sq is the Dirac mass at point (pay attention that our definition of the Green 
function obeys the convention used in probability theory: it is thus reversed in time 
in comparison with the standard notation used in the PDE literature). Following 
|H1 Th. 5, Sec. 5, Chap. 9], for a given (s,x) G [0, T) x (— oo, 1), the function 
(s, T] x (— oo, 1] 3 (t,y) i — y G(s,x,t,y) is also known to be the Green function of 
the adjoint operator 

d t -+d y [(b(y)+ae'(t)y]-^ y ; 
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with a Dirichlet boundary condition on [0, T] x {1}. In particular, G(s,x,t, 1) = 
and G(s, x, t, y) —> 5o(y — x) as t \ s. When xo = xo, we then set 

p(t,y) = G(p,x ,t,y), t e (0,T\, j/e(-oo,l]. (7.1) 

By [9] Th. 1.10, Chap. VI] (applied to the adjoint operator), we know that p(t,y) 
(as the solution to (3.4) under the current smoothness assumptions with Xq = xq) 
decays exponentially fast as t tends to and y stays away from Xq. This proves that 
p is continuous on any compact subset of ([0, T] x (— oo, 1]) \ {(0,Xo)}. 

Second Step. Still in the smooth framework, we now make the connection with the 
diffusion process x when X = Xq. Precisely, the point is to prove that the definition 
of p in (3.3) is coherent with that in (7.1). Put it differently, we must check that 
the right-hand side of (3.3) coincides with the definition of p in (7.1) when Xo = x o- 
Given a smooth function <ft : [0, T] x (— oo, 1] — > K, with a compact support, the 
analysis of the Green function in [9j Th. 1.10, Chap. VI] says that the PDE 

d s u(s, x) + (b(x) + ae'(s))d x u(s, x) + ^dl x u(s, x) + <p(s, x) = 0, 

for (s,x) G (0, T] x (— oo, 1), with u(T,x) = 0, x G (— oo, 1), as initial condition and 
u(s,l) = 0, s G [0, T], as Dirichlet boundary condition, admits a (unique) classical 
solution 



u(s, x) 



\ \ G(s,x,t,y)4>(t,y)dydt, sg[0,T), x < 1. 

J s J — oo 



(7.2) 



Moreover, u is bounded and continuous on [0, T] x (— oo, 1] and is once continuously 
differentiable in time and twice differentiable in space on [0, T) x (— oo, 1). Therefore, 
we can expand (m(Mt 1; XmtJ)o < t < r by Ito's formula, with T\ = inf{t ^ : Xt ^ 1}- 
We then have the well-known representation formula: 



u(Q,x ) 



E 



<f)(t,xt)dt. 







By equalizing (7.2) and (7.3), we deduce that 

TAn pT pi 



E 



4>(t,xt)dt 



o 



o 



p{t,y)<f>{t,y)dydt. 



(7.3) 



(7.4) 



Writing 



E 



TAti 



(f)(t,xt)dt 



[ E[^t,Xt)l{t< Tl} }dt 
Jo 



T pi 



<j)(t,y)F(xt edy,t< n)dt, 



we deduce that 

«T pi 







(j)(t,y)F( X t edy,t< n )dt 



p(t,y)<j>(t,y)dydt, (7.5) 



so that (3.3) holds in the smooth setting. 
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In the same framework, we then prove (3.5). The cumulative distribution function 
of 7~i is given by 



P( T1 < t) = 1 - P(n > t) = 1 - / 

J — CK 



p(t,y)dy. 



By [9 S Th. 1.10, Chap. VI], we can differentiate the above expression with respect 
to t and exchange the derivative and the integral. From (3.4), we deduce: 



d 
dt 



P(n < t) 



d t p{t,y)dy 

oo 

d y ([ae'(t) + b(y)]p(t,y))dy 



1 r 1 

2 y d 2 yy p(t,y)dy, £ > 0. 



Again by |9j Th. 1.10, Chap. VI], we know that both p(t,y) and d y p(t,y) tend to 
exponentially fast as y — > — oo. So, using the boundary condition p(t, 1) = 0, we 
obtain ( |3.5 ). 

Third Step. We now aim at proving the same results, still with % — x o, but 
under the original assumptions on b and e. The strategy is to use a mollification 
argument. In order to do so, we must prove that, in the smooth setting, p and d y p 
can be bounded by constants that only depend upon T, A, K, A and Xq, where A 
is such that sup 4 ^ T |e'(i)| ^ A. 

We thus go back to the case when b and e are smooth and bounded with bounded 



derivatives of any order. By Proposition 3.7, for t G [0, T] small enough, we already 
have a Gaussian bound for d y p(t, y) in terms of T, A, K and A only (notice that 
the argument applies since we know that p indeed satisfies PDE (3.4) in the smooth 
setting). With the same notation as in the statement of Proposition 3.7 (using in 



addition Proposition 3.5), the bound is of the form 



c 

\d y p{t,y)\ < -exp 



t 



|g° - y\ 
ct 



(7.6) 



for t G (0,5], where (£t°)t>o i s given by (3.14) and C and 5 are constants that 
depend on T, A, K and A only. Plugging this bound into (3.10) and repeating the 



Gaussian convolution argument used in (3.31), we deduce that, for t G (0,6] 



C 

P(t, y) < -j= exp 



\C t °-y\ 
ct 



(7.7) 



up to a new value of C. 

Actually, (7.6) and ( 7.7[ ) can be seen as bounds for the Green function G and 
its derivative in small time since p{t,y) = G(0,xo,t,y). In the same way, we could 
prove similar bounds for G(s,x,t,y) and d y G(s,x,t,y) when t — s G (0,5]. By a 



standard chaining argument, we then deduce that (7.6) and (7.7) are valid on the 
whole [0, T], for a possibly new value of C. Indeed, for a given t G (0,T], we can 
consider a sequence = to < t\ < ■ ■ ■ < t^ = t such that U+i — ti ^ 5. Then, by 
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using the Markov structure of x, we have: 



p(t,y) 



f* f N 

/ / I I G(ti^x, ti, Zi-i, Zi)dz\ . . . cfejv-i, 

J0 Ji-oo,!}"- 1 ~_i 



with the convention zq = and zn = U- Noticing that iV can be assumed to be 
bounded from above and using again a Gaussian convolution argument, we deduce 



that (7.6) and (7.7) can be extended to the whole (0, T\. 

Fourth Step. We still assume that the coefficients are smooth and bounded, with 
bounded derivatives of any order. By the third step, we are then able to reduce the 



PDE (3.4) to an heat PDE with a non-trivial source term: 



d t p(t,y) - l -d 2 yyP {t,y) = -d y [(b(y) + ae'(t))p(t,y)], t G (0,T], y < 1. 

For any compact subset JC C (0,T] x (— oo, 1], we can consider a smooth cut-off 
function 77 : [0, T] x R — > R + matching 1 on JC and vanishing outside another 
compact subset JC 1 C (0,T] x R, K C JC' . Then, the function (0,T] x (— 00, 1] 3 
(t,y) (—7- [i]p](t,y) satisfies the heat equation: 

(ft-^)M(*.v) = /(*,v). ( 7 - 8 ) 

with [?7p](0, •) = as initial condition and with [r)p](t, 1) = and [i]p](t, y) = for 
t G [0, T] and \y\ large enough as boundary conditions, where / is a smooth function 
which can be bounded in terms of T, A, rj, K and A only. By (HI Th. 4, Sec. 2, 
Chap. 7], we deduce that t]p and d y [qp\ are (1/4, l/2)-H6lder continuous in (t, x), 
the Holder constant depending on T, A, 77, K and A only. Therefore, on any compact 
subset K of (0,T] x (—00, 1], p and d y p are (1/4, l/2)-H6lder continuous in (t,x), the 
Holder constant depending on T, A, /C, K and A only. By the same argument, on 
any compact subset /C of [0, T] x X, X being a finite interval included in (—00, 1] not 
containing xq, p and d y p are also (1/4, l/2)-H6lder continuous in (t, x), the Holder 
constant depending on T, A, /C, K and A only. Indeed, in such a case, p(t, y) and 
d y p(t, y) tend to as t tends to and y stays in X, so that the compact support /C of 
r) can be assumed to be included in [0, T] xR (and not necessarily in (0, T] xl). In 
the end, we deduce that, on any compact subset JC of ([0,T] x (—00, 1]) \ {(0,Xq)}, 
p and d y p are (1/4, l/2)-H6lder continuous in (t,x), the Holder constant depending 
on T, A, JC, K and A only. 

In order to tackle the second-order derivatives in space, we assume that fC' C 
(0, T] x (—00, 1), KJ being the support of 77. Then, the function (0, T]xR9 (t, y) h-> 



[i]p](t,y) (with [rip\(t, y) = for y > 1) satisfies (7.8) on the whole (0,T] x R, so 
that it can be represented as a standard Gaussian convolution. Then, by Calderon 
and Zygmund estimates, see [TT?| Eq. (0.4), App. A], for any <; ^ 1, the L^([0,T] x 
R, dt (g) dy)-noims of d t [i]p\ and c^^p] are bounded in terms of A, rj, K, A and T 
only. Therefore, on any compact subset K. of (0,T] x (—00, 1), for any q ^ 1, the 
L ? (/C,<it (g) <i?/)-norms of c^p and <9 ro p are bounded in terms of T, A, JC, K and A 



only. For example, when JC is a cylinder of the form [S,T] x [y — l,y+ 1], (7.6) and 
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a 



f,<5 



ds, 



(7.9) 



(7.7) say that 

/ (\d t p(t,z)\" + \^ yy p(t,z)\")dtdz^C, >s fexp 

for a constant C fj< 5 that is independent of xq and y. 

Fifth Step. We now have all the required ingredients to go back to the original 
framework. The point is to approximate b and e by two sequences (b n ) n j> i and 
(e n ) n > i (for the topology of uniform convergence on compact sets) that satisfy the 
previous smoothness conditions. The associated solutions to the PDE ( 3.4[ ) are 
denoted by (p n ) n ^i- On any compact subset K C ([0, T] x (— oo, 1]) \ {(0,x )}, 
the sequences (p n ) n ^ i and (d y p n ) n i are uniformly bounded and (1/4, l/2)-H6lder 
continuous in (t, x). Therefore, there exists a subsequence (<£>(n)) n j>i such that 

^ i are uniformly convergent on compact subsets of ([0, T] X 
(— oo, 1]) \ {(0, Xq)}. Similarly, we can assume that the sequences (dtlp^^n ^ i and 
(dyy\P ))n > i are weakly convergent for the L ? (/C, dt <g> dy) topology, for any ^ ^ 1 
and any compact subset /C C (0,T] x (— oo,l). The limit function of the sequence 
(pv(™)) n ^ 1 i s denoted by p: clearly, it is a solution of (3.4) (in the Sobolev sense), 
with p(t, 1) = for t > 0, as Dirichlet boundary condition. By (7.7), it tends to 



as t tends to and y stays away from xq. Moreover, d y p exists and is continuous on 
([0,T] x (-oo,l])\{(0,xo)}. 

n ^ i converges toward 



Then, we can prove (3.3). Indeed, we know that (p 



p uniformly on compact subsets of ([0, T] x (— oo, 1]) \ {(0,a;o)}. Moreover, it is 
standard to prove that the solution x n to the SDE ( |3.2 ), but driven by b n and e n 
instead of b and e, converges in law (as n — > oo) toward \ on the space C([0,T],R) 
of continuous functions on [0, T] endowed with the uniform topology. Denoting by 
rf := inf{t ^ : Xt ^ 1}) we claim that the pair (r™ A T,x n ) converges in law 
toward (ti A T, x) on the product space K x C([0, T], R) endowed with the product 
topology. Indeed, the mapping 

C([0,T],R) 9i4 (inf{t ^ : x(t) ^ 1}) AT 

is continuous at any path x satisfying the "crossing property": 

(3t E [0,T) : x(t) ^ 1) (Vp ^ 1, 3t p E (t,t + l/p) : x{t p ) > l), 



(7.10) 



and ( 7. 10 ) holds for a.e. trajectory of x because of the Brownian part in the dynamics 



of X; see (3.2). Recalling that (7.4) holds true in the smooth setting, that is with 



p and Ti therein replaced by p n and r™, n ^ 1, we can pass to the limit along the 
subsequence (yj(n)) n ^ i (to pass to the limit in the right-hand side of ( |7.4 ), use the 
boundedness of the support of and the uniform convergence 

of (p&( n ») n>1 toward 

p on compact subsets of (0, T] x (— oo, 1]). We deduce that (7.4) holds in the general 
setting as well, from which we deduce that (7.5) also holds in the general setting, 
provided Xo 



X . 



By the same approximation argument, we can prove that (3.5) holds in the general 
setting as well, provided Xo 



x . 



Below, we shall denote p x ° for p to indicate the dependence upon xq. By (3.3), 
p x ° depends upon xq in a measurable way. 
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Sixth Step. To complete the proof, it remains to discuss what happens when xo 
does not reduce to a Dirac mass. The point is then prove that 



p(t,y) 



p x (t,y)F( Xo edx), t^O, y<:i, 



(7.11) 



is the right candidate for solving the Fokker-Planck equation and for making the 
connection with \. 



We first check that the definition (7.11) makes sense when t > 0. By applying 



the Markov property for x, this will directly prove (3.3). From (7.7), we know that 



p x (t,y) < (C/t 1/2 ) exp[-|£f - y\ 2 /(Ct)]. From ( [3.19P , we can find some C > such 
that 



cxp 



\tf-y\ 



^ I exp 

' — oo 



ct 

\*-e t \ 2 



p(xo e dx) 



C't 



= exp 
^ exp 



C't 

\i-e t \ 2 

C't 



P(xo e cfe) 
1 2(x-^t) 



C7't 



exp 



\x-e t \ 2 

C't 



(7.12) 



P(xo < x\dx 



■ i 2\ X - - xo\) 

1 , ^(2^) CXP 



is-&i 2 



which is finite since E|xo| < 00 > the third line following from an integration by parts 
and the last one from Markov's inequality. This also proves that p(t, y) — > when 
y — > — oo. 



The bound (7.12) can be easily modified into a domination argument for p x (t, y). 



for (t,y) in compact subsets of (0,T] x (— oo, 1]. By the Lebesgue Dominated Con- 
vergence Theorem, this proves that p is continuous on (0,T] x (— oo, 1]. When 
Xo ^ 1 — the same domination argument shows that continuity holds on any com- 



pact subset of ([0,T] x (— oo, 1]) \ ({0} x (— oo, 1 — e]), since the integrals in (7.12) 
run in that case from — oo to 1 — e. 

Using ( 7.6[ ), we can prove in a similar way that p is continuously differentiable in 
y on (0, T] x (— oo, 1], with 



d y p(t, y) 



I d y p x {t,y)F( X o e dx), t > 0, y < 1, (7.13) 

J — oo 



and that d y p(t, y) — > when y — > — oo When xo ^ 1 — continuous differentiability 
holds on any compact subset of ([0,T] x (— oo, 1]) \ ({0} x (— oo, 1 — e]). 

By combining the same domination argument with (7.9), we can also prove that 
p admits Sobolev derivatives of order 1 in t and of order 2 in y in any L ? , c ^ 1, on 
any compact subsets of (0, T] x (— oo, 1), the derivatives being given as the integrals 
of the derivatives of p x with respect to the law of Xo- It is then plain to check that 



p satisfies (3.4) 
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The last point is thus to check that p satisfies (3.5). By Markov's property, 
P(Ti<t)= I F(r 1 ^t\xo = x)F(xoedx) 

J — oo 



d y p x (s, l)ds 



(7.14) 



P( Xo e dx). 



In order to prove (3.5), we must exchange the two integrals. To exchange the two 
integrals, we specify the upper bound for \d y p x (t,l)\. Using (7.6) and (3.19) and 
recalling that = 1, we deduce that \d y p x (t, 1)| ^ (C/t) exp[— 11 — x\ 2 /(Ct)], from 
which it holds that 



1 r t 



oo JO 



\d y p x (s,l)\dsF(xo G dx) 



oo JO 



C 



exp 



+ 



1 rt 



l-e JO 



c 



exp 



[1-xj 2 
Cs 

|1 — x\ 

Cs~ 



ds P(xo e dx) 



dsF(xo e dx). 



(7.15) 



The first integral in the right-hand side is clearly bounded. To tackle the second 
one, we make use of the assumption P(xo £ dx) ^ /3(1 — x)dx for x G (1 — e, 1]. We 
obtain 



l-e Jo 



C 



exp 



x 



Cs 



ds¥(xo edx) ^ /3 



6 



^0 



Co: 



exp 



x 



Cs 



ds dx < oo. 
(7.16) 



Going back to (7.14), we obtain 

• i 



P (rj < *) 



1 



o 



d y p x (s, l)P(xo G ofo^ofs 



d y p(s, l)ds, 



the last part following from (7.13). Since the mapping (0, T] 9 t i— )• d y p(t,l) is 
continuous, we deduce that the mapping (0, T] 9 t \-¥ P (n ^ t) is continuously 
different iable. It then remains to check the continuous differentiability at time t = 0. 
To this end, it is sufficient to prove that d y p(-, 1) is continuous at t = 0. Following 
(7.15) and (7.16), it is clear that, for any e' G (0,e), 

limd„p(t, 1) = lim / d y p x (t, 1)P(y G rfx) = lim / ^(t, lWaOdz, (7.17) 

i\0 t\0 < / 1 _ e , i\0 y i _ e , 

where poO 3 ?) — {d/dx)W{xQ £ rf^)) which is assumed to make sense on (1 — e', 1). By 
Proposition |3.7[ we know that, for some constant C > 0, 
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which proves, by using the bound po(x) ^ 0(1 — x) and by letting e' tend to in 
(17.171), that 



limcLp(t, 1) = lim lim / d v q(t, x,l)p { 

t\0 y e'\0t\0 J x _ e , J 



x)dx. 



(7.18) 



We then introduce the killed Gaussian kernel with reflection at point 1: 



q{t,x,y) 



lixt 



exp 



\x - y\ 



exp 



\x-{2 



2t J V 2t 

Now a formula like (3.10) for expressing d y p^(t,y) in terms of d y q(t,x,y) can be 
proven for expressing d y q(t,x,y) in terms of d y q(t,x,y) (replacing ae'(s) + 6(1) by 
b(z) — b(l) in ( 3.10[ )). Recalling \b(z) — b(l)\ ^ K\z — 1|, we obtain (for a possibly 
new value of C which is allowed to increase from line to line) 

\d y q(t, x, 1) - d y q(t,x, 1)| 



^ C 
< C 

^ c. 



-oo 
t /■! 



1 ^(1 ~ * ) 

s(t - s) 
C 



exp 



( 


2 ^ 


exp 




1 - 


3 









S (t-s)V2 



exp 



dzds 



Cs 



dzds 



Recalling (7.18), we deduce that: 



limcLx>(t, 1) = lim lim / d v q(t, x, l)pn{ 

t\0 E'\0t\0,/ w 



x)(ia;, 



p (x)dx, 



= - lim lim 2r 1 / 2 / fi 1 ' ( "~T7T 

where ^ stands for the standard Gaussian kernel. Since p is assumed to be differen- 
tiable at 1, with po(l) = 0, we can write po(x) = (x — l)p' (l) + o(x — 1), where o(-) 
stands for the Landau notation. We deduce that the limit of d y p(t, 1) as t \ must 
be Po(l). Following (7.15) and (7.16), we have an explicit bound for the supremum 
norm of d y p(-, 1) in terms of a, f3, b, T and the supremum norms of e and e' only. 
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